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Semiclassical spin-coherent kinetic equations can be derived from quantum theory with many differ- 
ent approaches (Liouville equation based approaches, nonequilibrium Green's functions techniques, 
etc.). The collision integrals turn out to be formally different, but coincide in textbook examples as 
well as for systems where the spin-orbit coupling is only a small part of the kinetic energy like in 
related studies on the spin Hall effect. In Dirac cone physics (graphene, surface states of topological 
insulators like Bii-^Sb^, Bi2Te3 etc.), where this coupling constitutes the entire kinetic energy, the 
difference manifests itself in the precise value of the electron-hole coherence originated quantum 
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correction to the Drude conductivity ao ~ x^ F - ^he leading correction is derived analytically 
for single and multilayer graphene with general scalar impurities. The often neglected principal 
value terms in the collision integral are important. Neglecting them yields a leading correction of 
order (Ikp)' 1 , whereas including them can give a correction of order (£k-p)°. The latter opens up 
a counterintuitive scenario with finite electron-hole coherent effects at Fermi energies arbitrarily far 
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above the neutrality point regime, for example in the form of a shift 8a ~ %- that only depends on 
the dielectric constant. This residual conductivity, possibly related to the one observed in recent 
experiments, depends crucially on the approach and could offer a setting for experimentally singling 
out one of the candidates. Concerning the different formalisms we notice that the discrepancy be- 
tween a density matrix approach and a Green's function approach is removed if the Generalized 
Kadanoff-Baym Ansatz in the latter is replaced by an anti-ordered version. This issue of Ansatz 
may also be important for Boltzmann type treatments of graphene beyond linear response. 



I. INTRODUCTION 



Since the first isolation of graphene in 2004 pQ the electrical conductivity of this system has attracted huge attention. 
To a good approximation the electrons can be described as massless 2d Dirac electrons for which the spin-orbit in- 
teraction that yields the characteristic Dirac cone is given by the pseudospin derived from the bipartite honeycomb 
lattice. [51 13] The conical electron and hole bands touch at the Dirac points. The Brillouin zone contains two inequiv- 
alent degenerate Dirac points — K and K' — that give an additional valley index. Finally there is the ordinary electron 
spin. To first approximation, the conductive properties of graphene involve only the pseudospin in a nontrivial way. 

The main focus has been on undoped graphene, with the chemical potential exactly at the degenerate Dirac points. 
This regime of chemical potential close to zero — the Dirac regime — hosts the most exotic features, for example the 
finite conductivity minimum at seemingly zero charge carrier density. ^1 E] Quantum effects due to electron- hole 
coherence (that is, pseudospin coherence) like Zitterbewegung can dictate the observed conductivity even to lowest 
approximation.^ We refer to ref. [7j for a review on early work on the Dirac regime. 

Away from the Dirac regime, with a large enough charge carrier density there is a crossover into the Boltzmann 
regime Ik-p 3> 1 (with £ the mean free path and hk-p the Fermi momentum) . Here, the conductivity can be understood 
to lowest order without taking into account quantum effects such as electron-hole coherence and are therefore more 
intuitive. The crossover between the two regimes has recently been studied numerically. PHTU] 

In many experiments on graphene samples on substrates [U [5J 111H13) the dc conductivity in the Boltzmann regime 
is observed to be linear in the electron (hole) density (see fig. [IJ. This results in the characteristic V-shape in 
the conductivity as a function of gate voltage. The linearity is less pronounced in suspended graphene where the 

concentration of charged impurities in particular has been reduced. [141 115] Theoretically the linear behavior is well 
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described by the ordinary Drude conductivity 4<7o = 4 ^tk-p derived from a Boltzmann equation with four degenerate 
(valley and real spin) incoherent bands, provided that screened charged impurities are assumed. |16H18j Point-like 
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FIG. 1: Sketch of the dc conductivity a in monolayer graphene as a function of the electron density n as observed in some 
experiments [4l l5l 1111413] . in particular including a residual conductivity Sa as observed by Chen et al. [13j . At the neutrality 
point n = the charge density is zero and one would expect zero charge carrier density at low temperatures. One of the 
surprises of graphene is the conductivity minimum <r m ~ e 2 /h. In the Dirac regime n ~ (red/gray) the usual criterion 
IkF 3> 1 for a Boltzmann treatment is not satisfied (at least not for screened charged impurities) . Away from the Dirac regime 
a semiclassical approach should usually apply. The observed linear dependence in the Boltzmann regime (black) is described 
by the Drude conductivity ctq = 2e 2 Iky/h oc |n| for screened charged impurities. (Point-like impurities, in contrast, yield 
(To oc \n\°. The resistivity due to them starts to compete with that of screened charged impurities for |n| large enough.) Effects 
of electron-hole coherence do not enter into the Drude conductivity but reveal themselves only in quantum corrections of higher 
order in (Ik-p/ 1 ■ A contribution (^/cf) _1 can explain the initial convexity as one approaches the Dirac regime. A contribution 
(£fcp)° enters as a constant shift in the Boltzmann conductivity and contributes thus to the residual conductivity Sa. The 
latter is directly read off by linear extrapolation. A precise knowledge of other contributions to Sa (e.g. weak localization, . . . ), 
allowing precise estimate of the electron-hole coherent contribution, which in turn would single out one of the many possible 
derivations discussed in the present paper. With one except ion 1 10] the contribution to the conductivity minimum depends 
also in general on the approach. However, these differences can probably not be studied as cleanly as those in the residual 
conductivity. We want to stress that the residual conductivity is inherent in the Boltzmann regime whereas the minimum 
conductivity is inherent in the Dirac regime. 



impurities on the contrary lead to a Drude conductivity that is independent of charge carrier density and therefore 
this model fails even on a qualitative level. [19] As one expects the charge carrier to be mainly of one type (either 
electrons or holes), the pseudospin band index can be left aside. With interband-coherences neglected the collision 
integral contains only transition rates between energy eigenstates. Such rates are simply derived with Fermi's Golden 
rule. The specific properties of massless Dirac electrons enter merely into the transition rates as a spin-overlap 
factor due to the chirality of the eigenstates and as a Fermi momentum dependent (and therefore electron density 
dependent) Thomas-Fermi momentum j»tf due to the linearity of the spectrum. (In a 2DEG with quadratic dispersion 
the screening length h/p^-p is in contrast independent of the Fermi momentum.) 

In the Boltzmann regime £kp ^> 1 the electron-hole coherent features of Dirac electrons only manifest themselves 
if one goes beyond the lowest-order result <7o oc Ik? to look for quantum corrections. A Boltzmann approach to these 
quantum corrections requires a kinetic equation that is quantum coherent in band indices, for graphene the pseu- 
dospin index. [2"UH2"2"] Interband coherent collision integrals are beyond the application range of Fermi's Golden rule. 
To access the "transition rates" involving the interband components (a.k.a. off-diagonal, (pseudo)spin-precessing or 
Zitterbewegung components) in the collision integral one typically has to resort to a fully quantum coherent theory 
and then derive a Boltzmann equation by a semiclassical expansion in the space and time variables but not in the 
spin, which is to be kept quantum coherent. Common approaches for deriving quantum coherent kinetic equations 
are density matrix approaches and nonequilibrium Green's function approaches, both with further subdivisions. The 
former start with a single-time density-matrix-type state variable p{x\, t\) governed by a quantum Liouville equa- 
tion (von Neumann equation). The latter start with a double-time correlator G < (xi,ti, £2,^2) governed by dynamic 
equations derived from e.g. the Kadanoff-Baym equation or the Keldysh equation. At a later stage some approxi- 
mation has to be invoked to recover single-time equations. This choice of approximation is the problem of Ansatz 
which will be discussed later. All the mentioned treatments of pseudospin-coherence induced quantum corrections 
in graphene use the density matrix approach. In other contexts involving spin-coherent Boltzmann approaches to 
transport in the presence of spin-orbit coupling such as the spin Hall effect, Green's functions derivations have also 
been used. |23ff27] However, in these cases transport was addressed to leading order rather than quantum corrections 
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One of the two central questions of this paper is whether the two approaches are equivalent in general and in particular 
when addressing quantum corrections to the graphene conductivity. 

The general small expansion parameter for quantum corrections in a Boltzmann regime is (ikp)^ 1 . One example 
is the famous result by Gorkov et al. [28] for the correction Serf <jq = (ir£k-p)~ 1 log loqt due to weak localization. 
The quantum correction to the graphene Drude conductivity due to electron-hole coherence is no exception. [61j The 
second central question addressed in this work is whether the series of quantum corrections in powers of {£kp) _1 
starts at the order (£k F )~ 2 a or already the order (£kp)~ 1 a . A correction 6a ~ ^(ikp) -1 to the Drude conductivity 
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~ TT^f would depend on impurity concentration and impurity strength. For screened charged impurities in 
graphene it would increase with decreasing electron density n, which would explain the onset of convex behavior 
of the conductivity as one approaches the Dirac regime. Far away from the Dirac regime the contribution becomes 
arbitrarily small. A contribution of the lower order 6a ~ tt{^f) would be independent of the impurity density and 
impurity strength. At least in monolayer graphene it would also be independent of the electron density, which would 
appear as a constant shift of the Drude conductivity as illustrated in fig. [I] Electron-hole coherent effects would then 
remain finite arbitrarily far away from the Dirac regime, which appears rather counterintuitive. Refs. [201 121| found 
the leading correction to the dc conductivity to be of the order (ikp)^ 1 within a pseudospin-coherent Boltzmann 
approach. (62) Trushin and Schliemann|21) found a leading correction of the same order |63j although their approach 
was qualitatively different in that they, in contrast to Auslender et al., discarded the principal value terms in the 
pseudospin coherent collsion integral. (See also ref. [21].) Both papers found the result to be ultraviolet divergent for 
point-like impurities. |64| Culcer and Winkler[22 studied screened charged impurities. They too neglected principal 
value terms. They solved their spin-coherent Boltzmann equation only up to order (£kp)°, wherefore the previously 
found corrections ~ (^fcp) -1 were out of reach. However, solving the Boltzmann equation by introducing a frequency 
dependence lead to a quantum correction 6a(u)) ~ e 2 /h for non-zero temperatures. As this correction of order (£kp)° 
vanishes in the zero-frequency limit their result is not inconsistent with that of Trushin et al. or, at first sight, that 
of Auslender et al. . 

An important contribution of the paper of Auslender et al. |20j was that they realized the importance of principal 
value terms and included them in their collision integral for graphene. However, they did not discuss the physical 
origin or necessity of these terms that so severely complicated their analysis. The subsequent papers jHHIH] neglected 
them without much comment and this seems to be the rule also in other contexts where Boltzmann approaches are 
applied to spin-orbit coupled systems, as for example in the spin Hall effect . [2"3Tf2"7] In other fields the meaning of these 
terms in a Boltzmann treatment is rather clear (for the context of strongly interacting fermions see e.g. Lipavsky et 
al. [30]). As this does not appear to be the case in the spin-orbit context, we will in the technical sections discuss the 
physical origin of these terms and thus motivate why we insist on keeping them although they make the interpretation 
of the Boltzmann equation more elusive and also considerable increase the technical challenge of solving these integro- 
differential equation. A witness of this increased complexity is the mathematical tour de force of Auslender et al. |20j , 
and that treatment still applying only to the simplifying case of point-like impurities. 

In trying to answer the two mentioned central questions the present paper has two aims. One aim is to (re)derive 
a pseudospin-coherent Boltzmann equation and to give a systematic derivation of the leading quantum correction 
by including the principal value parts, but nonetheless consider screened charged impurities which are relevant to 
graphene. We can consider a simpler problem than that of Auslender et al. by using an iterative scheme in the 
spirit of Culcer et al. rather than attempting to solve the equations to all orders in (^fcp) -1 in one go. To illustrate 
the importance of the principal value terms we also solve the same Boltzmann equations with principal value terms 
neglected. This problem, in contrast, is also easily solved non-recursively. [27] Since we formulate the Boltzmann 
equation for a spin-orbit coupling of arbitrary winding number N the results also apply to certain models [31-34 for 
bilayer and multilayer graphene. Quantum corrections in the bilayer case N = 2 have been considered by Culcer et 
al. [35] with a similar approach as in ref. [52]. We will refer to N = ±1 as the "monolayer case" and to \N\ ^ 1 as 
the "multilayer case" . (Actual iV-layer graphene hamiltonians can, however, depending on the stacking, be written 
as tensor products of lower-JV hamiltonians, in some cases including the monolayer N — ±1 hamiltonian433, 34 ) Wc 
find the following results: 

1. With principal value terms neglected, the first quantum correction is of order (^f) -1 both in the dc conductivity 
as well as in the dissipative ac conductivity. Higher order corrections beyond the order (^fcp) -1 are absent in 
the monolayer case \N\ = 1, but present to infinite order in the multilayer case \N\ ^ 1. 

2. With principal value terms included, the leading quantum correction to the dc conductivity is of order (£kp)°. 
In particular it is independent of the impurity concentration and of the impurity strength. For screened charged 
impurities at a negligible distance from the graphene plane, the correction depends only on the dimcnsionless 
parameter fcp/^TF (with hk-p the Fermi momentum and 1/fcxF the Thomas- Fermi screening length). In the 
monolayer case TV = ±1 (with &tf &f) the correction is also independent of kp, that is, of the electron 
density, and depends only on natural constants and the dielectric constant. 
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3. For screened charged impurities the corrections are convergent. For point-like impurities the corrections are 
convergent in the multilayer case \N\ > 2 but require an ultraviolet cut-off in the monolayer case N — ±1. 

4. For point-like impurities in the multilayer case \N\ ^ 1 the contribution from principal value terms vanishes 
trivially to orders (£fc F )° and (ffcp) -1 - The leading quantum correction is then the correction of order (Ikp)^ 1 
derived with principal value terms neglected. 

The first conclusion is in qualitative agreement with Trushin et al. [2T]. The ac result does not correspond to the 
frequency-dependent correction ~ (ikp) found by Culcer et al. [33] [35]. The second conclusion appears to disagree 
with Auslender et al. [20] who present a leading quantum correction of order (Ihp) -1 .\Q5\ The third conclusion is in 
agreement with all the previous studies. The fourth conclusion is consistent with the recent paper ref. |10j.|66j We 
wish to stress that our electron-density-independent contribution ~ (£k-p)° to the conductivity has a different origin 
than those discussed in previous studies [2lJ [22] . 

Electron-hole coherent quantum correction gives with the shift (£kp)° one mechanism for generating a residual 
conductivity(see fig.[l]) . Our shift could therefore play a part in the residual conductivity observed in the experiments 
by Chen et al. [12] on monolayer graphene. This residual conductivity — estimated to be 2.6e 2 //i — is observed to be 
surprisingly constant, depending at the most weakly on impurity density, in contrast to the conductivity minimum. 
Our contribution to the residual conductivity is independent of the impurity density. However, can in principle 
depend on the dielectric environment of the graphene sample. There are experiments varying the dielectric constant 
for example by coating the sample with ice|36j. However, here the extraction of the residual conductivity appears to 
be more precarious, and it is too early to say if there is any relation with the dielectric behavior of our contribution. 

The second aim of the paper is to compare different derivations of collision integrals for graphene and to extend 
the repertoire to Green's functions approaches. Double-time approaches require, however, the choice of an Ansatz. 
We will see that different choices lead to qualitatively different general collision integrals. Some of our results are: 

1. The standard choices of the Kadanoff-Baym Ansatz (KB A) [37] or the Generalized Kadanoff-Baym Ansatz 
(GKBA) [38] do not give the same general collision integral as the single-time density matrix approaches. 

2. We propose an alternative Ansatz — AA ( "anti-ordered Ansatz") — for which the translation between density 
matrix and Green's functions approaches can be established. 

3. For the symmetrized KB Ansatz (SKBA) given by the sum of GKBA and AA the principal value terms can be 
completely absent. 

4. For spinless electrons and to zeroth order in gradient expansion the general collision integrals of all approaches 
coincide. This applies even in the presence of spin-orbit coupling provided the latter is small. The difference 
between collision integrals enters into the principal value terms, which vanish for zero spin-orbit coupling (to 
zeroth order in gradient expansion). The difference also enters into the delta function terms, but only in parts 
that capture second and higher order effects in the spin-orbit interaction. 

5. Both with and without principal value terms, the leading quantum correction to the graphene Drude conductivity 
depends generally on the approach. The difference can be quantitative (for example a factor 3 for point-like 
impurities in the monolayer case \N\ = 1 with principal value terms neglected when comparing the result derived 
with a the density matrix approach to the result derived with a Green's function approach implementing the 
GKBA) or, when principal value terms are included, even qualitative (of different order in (ik-p)" 1 or of opposite 
sign, in the case of point-like impurities in monolayers \N\ = 1 of different ultraviolet divergent behavior). 

The exception to the last conclusion is the case of point-like impurities in a bilayer | N\ — 2. In this case the principal 
value terms vanish in all approaches and the remaining delta function terms of the collision integral coincide. Therefore 
the density matrix approach of Ref. [TU] is not concerned by the ambiguities unravelled in the present paper. 

Conclusion 4, partly also discussed in ref. [27] , might explain why the difference between a density matrix approach 
and a Green's function approach with the the GKB Ansatz does not seem to have been an issue of debate [57], even 
in similar contexts like spin-coherent Boltzmann treatments of the spin Hall effect, where there is a bigger variety 
of applied approaches, see refs. [2"3Tt27| . With the focus on spin and charge currents to leading order, the typical 
neglect of principal value terms seems acceptable. Furthermore, the energy b of the spin-orbit coupling is typically 
assumed to be small compared to the spin-diagonal part e ~ p 2 /2m of the kinetic energy to motivate an expansion 
in the collision integrals to first order in bp/e^p, see [HI 123 127] ■ The difficulty of solving the Boltzmann equation 
analytically increases rapidly with higher orders in bp/e^p an d appears to be beyond reach without an expansion. 
However, in the special case of eo = 0, the case of graphene, many simplifications occur which enable an analytic 
solution even for a strong spin-orbit coupling. The circumstance eo = is what allows us to analytically explore the 
differences in predictions between the different approaches. 
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Apart from the issue of the principal value terms the common framework of approximations of the related works 
are adopted in the present paper. In particular we assume a low impurity concentration, restrict ourselves to the 
lowest order in Born approximation and neglect weak (anti-)localization corrections. This framework also included 
two further approximations: 

1. Terms involving the interactions (self-energy terms) are taken only to zeroth order in gradient expansion while 
the free part of the kinetic equation is expanded to first order in gradients (in order to recover the driving 
terms). 

2. The Wigner transformed density matrix of the electron p(p,x,t) ~ / dui G < (p, u, x, t) is identified with the 
Landau quasiparticle distribution function f{p,x,t). The former is the quantity for which the kinetic equations 
are formally derived and in terms of which the currents are expressed, the latter the quantity which in equilibrium 
is given by a Fermi-Dirac distribution. 

Although each one of these approximations seems to be inconsistent at first sight, both are standard. They are 
perfectly fine for the derivation of a lowest order result like the Drude conductivity or the spin Hall current. These 
approximations break down as soon as one has to consider quantum corrections due to strong interactions or strong 
fields beyond linear response. (See refs. [301131].) However, it would be far too ambitious at present to also properly 
account for gradient corrections and the difference between p and / in our study. Also, such a treatment would not 
directly be concerned with the two mentioned purposes of our study, but rather with the global goal, namely that 
of a fully consistent and systematic Boltzmann treatment of the first quantum correction. This problem is work in 
progress that we hope to return to in the future. 

In this paper we have mainly graphene in mind. However, the results apply to any setting in 2d where the electrons 
are described by one or several independent Dirac cones. Our results should therefore also be relevant for some 
3d topological insulators like Bii-^Sba;, Bi 2 Te3, Sb2Te3 and Bi 2 Se3 with the 2d surface states described by Dirac 
electrons [4"Ull45| . see also ref. [46] . In particular, the occurrence of a single Dirac cone removes the problem of 
intervalley scattering. Thus, these systems should be a better application of those of our calculations that deal with 
point-like impurities, for which we would expect our assumption of negligible intervalley scattering to be invalid. 

The outline of the paper is the following. In sec. [H] we present the Wigner transformed hamiltonians that we 



examine in this paper. In sec. Ill the semiclassical distribution function and the Boltzmann equation for a spin-orbit 
coupled system are introduced. Derivations of collision integrals with different approaches are presented in sec. |IV| In 
sec. [V] the different general collision integrals are compared and simplified for the case of spin-orbit coupled systems, 
in particular graphene. Sec. |VI| deals with the solution of the Boltzmann equation neglecting principal value terms 
and the resulting dc conductivity. Sec. VII contains the ac current generalization thereof. In sec. |VIII we solve the 



Boltzmann equation including the principal value terms and present the resulting conductivity, both dc and ac. In 
the beginning of that section we also discuss the physical background of the principal value terms. At the end of that 
section we discuss whether finite electron-hole coherent effects far away from the Dirac point make sense. Conclusions 



and an outlook are given in sec. IX 



II. THE MODEL 



For a semiclassical Boltzmann description (see e.g. [39 ll47ti49] one typically uses the Wigner transformed one-particle 
Hamiltonian. For the spin-orbit coupled systems that we set out to study the hamiltonian in absence of impurities 
reads 

H(x,p,t) = e (k) + tr- b(k) + e<f>{x,t) (1) 

with e < and k(x,p,t) = p — eA(x,t). The energy bands in the absence of electromagnetic fields are e| — eq + sb 
with s = ± giving the sign of the spin along the spin quantization axis b, i.e. er ■ b\bs) = s\bs). We want to describe 
a 2d system with x and p chosen to lie in the x, y-plane. We set H = 1. 

The spinless part of the dispersion is given by eo. In the spin-orbit coupled systems studied in the intrinsic spin 
Hall effect, this is typically the quadratic dispersion eo = k 2 /2m* and usually consitutes the bigger part of the kinetic 
energy. For monolayer graphene, the spin-orbit coupling entirely constitutes the kinetic part, i.e. eo = 0. The spin is 
here the pseudospin given by the bipartite lattice. Additionally, the electrons have a valley index corresponding to 
the two Dirac cones, K and K', as well as the real spin index. The real spin index will be treated as trivial in the 
present paper. Furthermore, we will neglect inter-valley scattering to be able to treat each Dirac cone independently 
and therefore deal only with the 2x2 hamiltonian ([I]). This approximation should be fine if the disorder is smooth 
but should brake down if the impurity potential is too short-ranged, like in the extreme case of point-like impurities. 
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For the Dirac point K the spin-orbit coupling for the pseudospin is given by b = vpk (i.e. b = vpk and b = k) with 
tip ?s c/300 being the Fermi velocity. We consider the more general isotropic spin-orbit coupling b ~ b(k)b(0) with 
a winding number N given by b x + ib v = e l8o+lNe (80 is a constant). This includes the Dirac cone K' (b — vpk and 
b = (cos#, — sin#), i.e. N = —1) of monolayer graphene and the hamiltonian 

_ If (k x Tik y )i\ k 2 f e=F«'\ 
M ~ 2m* \(k x ±ik v f ) ~ 2m*\e ±aB ) [A) 

(b = k 2 /2m* and N = 2) studied in the context of bilayer graphene as well as similar hamiltonians studied in 
multilayer graphene. [3TH34] 

The total hamiltonian H to t = H + V includes an impurity potential V(x) — ^2 n u(x — x n ) of non-magnetic 
impurities at positions x n eventually to be averaged over. We distinguish between point-like impurities u kk > — const, 
and screened charged impurities in 2d with 

e 2 /«d ,„s 
Ukk ' = |fc-fc'| + fc TF (3) 

with the Thomas-Fermi momentum fcxF = ^^D(e-p) giving the range L s — 1/fcxF of the screened potential, is the 
dielectric constant. For a recent review on screening in graphene with an enlightening comparison of monolayers with 
bilayers and 2DEGs, see ref. [SO]. Here we will recollect a few facts that will be important for our later discussions. 

It is convenient to introduce the dimcnsionlcss parameter q s :— /ctf/^f characterizing the strength of the screening. 
An unscreened Coulomb interaction corresponds to q s = 0. In the opposite limit q s —¥ 00 the potential becomes 
almost angularly independent and behaves in some respects as a point-like impurity, although not in all aspects. In 
the monolayer case, for example, one has the unique situation that Dp cx hp, hence q s is independent of kp , implying 
that for short and long screening lengths alike does the potential behave as an unscreened Coulomb potential in 
that T t ~ 1 (fcp) oc Dpkp 2 cx fcj 1 and therefore cto ~ T tr (kp)ep oc fc| oc \n\. Thus, not even at strong screening does 
the situation turn into that of point-like impurities, where r^ 1 cx fep predicts a Drude conductivity independent of 
density. For graphene on Si02 substrates, the standard value is 

q s w 3.2 . (4) 

This suggests that screening is important (q s > 1) and when discussing the quantum correction (£kp)° we will argue 
that this correction for the screened potential with the given q s is nonetheless closer related to the results for point-like 
impurities than to those of an unscreened Coulomb potential. 

In bilayers and multilayers the situation is very different. Here q s decreases with fep and consequently with the 
density, just like in an ordinary 2DEG, but opposite to a 3DEG. Thus, the further we move away from the Dirac 
point, the weaker the screening and the stronger the effect of the interaction. Thus, when discussing the correction 
(^fcp) we expect that in the vicinity of the Dirac point the screened potential has more in common with point-like 
impurities where as far away from the Dirac point the potential has more in common with an unscreened Coulomb 
potential. 



III. SEMICLASSICAL DESCRIPTION OF A SPIN-ORBIT COUPLED SYSTEM 



In a Boltzmann picture the spatial degrees of freedom can be treated semiclassically. The treatment of the peudospin 
must, on the contrary, remain quantum mechanical in order to capture electron-hole coherence effects. The state of 
the system is described by the Wigner transform p a(J r(x,p, t) of the time-diagonal density matrix 

Paa'{xi,x 2l h) = (•^l,(x 2 ,t 1 ) 1 ii a (x 1 ,t 1 )) = G^ a , (x 1 ,t 1 , x 2 ,t 2 )t 2 =t 1 (5) 

for electrons with spin indices a =tz,iz- (See e.g. refs. [3T)I IT7H4T)] .) In the absence of interactions one can derive a 
Boltzmann equation for p by applying the Heisenberg equation of motion on ^(ii), then identifying t 2 = ti, Wigner 
transforming the result, and gradient expanding to first order. The approximation to stop at first order in gradient 
expansion is the semiclassical approximation, which assumes that the external perturbations, such as electromagnetic 
potentials, change negligibly on length and time scales of the de Broglie wavelength Ab and the time tb = Ab/uf- Upon 
Wigner transformation X(xi,ti,X2,t2) — > X(x,p,t,u>) (for the time-diagonal approaches the Wigner transformation 
only involves the spatial part, not time) the convolution product is translated into XY — > Xe^Y where the matrix 
product in spin remains but the convolution of time and real space variables is replaced by a Moyal product with a 
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Poisson-bracket-like gradient T>. In a gauge invariant treatment, valid when the electromagnetic fields are weak and 
vary slowly (see e.g. ref. |47|) . one introduces k(p, x,t) = p — eA and Q(uj, x, t) = ui — e<\> and lets {x, k, t, ui} become 
the new set of independent variables (i.e. d Xi k — 0). This changes the gradient into 

V = d Xi d ki - d k ~d Xi + du~o t - d t OQ + Ei(d Q ~[) ki - dk~Ou) + ^iji B % 9 k s ~o k, (6) 

with X*dY := (dX)Y and X~$Y := X(dY). 

From the matrix elements of the distribution function p, one extracts the densities and current densities of charge 
and spin. The matrix elements are conveniently expressed in the decomposition p = p + <x ■ p where er = (cr x , a y ,a z ) 

is the vector of Pauli matrices. Furthermore, we find it convenient to decompose the vector p = p^b + p^c + p z z in 
its components along the basis vectors b(9), z and c(9) — z x b(9) with deb = Nc , analagous to the cylindrical basis 
vectors k(9) := k/k, z and 0(9) := z x k(9) with dgk — 6. In the spin basis {| tz)i I iz)} one has 

p = p l + p h b • a + P£ c • a + Pz z • a= ( (p /^> ^ ) • (7) 

The charge density en and current density ej in phase space are derived from en — Tr (pdH/d(jj) and ej = 
-Tv(pdH/dA), which yields 

n(x,k,t) := Trp = 2p = n + +n~, 

ji(x, k, t) := Tr ( Vi />) = 2p d fci e + 2p • d fei 6 - n+v+ + n~vr + ^p £ 9 t ( > 

with i = x,y. Here we introduced the velocity matrices Vi := d ki H = dk^o + c ■ d^b. The spin-independent part of 
the velocity is dk£o =■ t>o- The band velocities are v s := dk£ s = (bs\v\bs) — v s k. The m£ra-band elements 



± 



= (b±\p\b±)=p ± Pi (9) 



give the density in each spin band s = ±. The mier-band elements (6 ± |p|b=p) = p z ± ipg are important for the 
coherent treatment of spin. These are the elements that oscillate in the occurrence of spin-precession. In the case 
of spin-orbit coupling the imaginary component pa appears in the last term of the current Q |68j . For the density 
matrix of a single electron this term would contain the oscillatory Zitterbewegung motion of the free spin-orbit coupled 
electron. In the statistical description given by the distribution function this oscillatory motion of the free particle 
states average to zero over time and is therefore absent in the equilibrium distribution function (p? q = p z q = 0). 

Throughout the paper we write the Boltzmann equation and the distribution function in the momentum-independent 
spin basis {| t*)>| -J-)}- Some studies (e.g. refs. [ST]) prefer to transform the Boltzmann equation into the 
momentum-dependent eigenbasis {|i>+), \b— }}. For example, the velocity matrix for graphene (eo = 0) 



Nb ( kd k b iO 



Nb 



cli 



v = a t (kd k + \0d )bi = ka- bd k b + 6 a ■ c— = ^ j (10) 

when written in the latter basis (the superscript "ch" indicating the chirality/helicity basis). For the distribution 
function the relation is 

(\ ch 
po + P 6 p, + m =pol + ^ _ pea o b + pz(7 ch 
Pz l Pc PO l J b / 

The Pauli matrices af 1 have non-zero momentum derivatives (dgcrf 1 = —Na y h , de<r y h = Naf 1 and cW£ h = 0). 
The eigenbasis has the attractive feature of allowing one to identify the intra-band densities n = P o ± pg as the 
diagonal matrix elements. However, this would not be of any advantage in our treatment of graphene, where we 
only need the spin components and need to distinguish between pg, pg and p z rather than to view the problem in 

terms of p_| (= pl + ), p++(= n + ) and p (= n~). The charge component po will in our problem be redundant and 

completely decoupled. 

The spin density, i.e. the polarization, is given by s M = ^Tr (cr^f) = p^ (with p = x,y,z). There is not a unique 
way to define the spin current because in presence of spin-orbit coupling, the real space spin polarization is not a 
conserved quantity. When band velocities coincide (v s = Vq) it is = p^vq. A common definition is 

3i = ^ Tr 0>{v;, p} = p^d ki eo + pod^fyu (12) 
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(with {A, B} = AB + BA). Notice that in the graphene case eo = one is confronted with the peculiar situation 
that the charge current is only determined by the spin density, whereas the spin current is only determined by the 
charge density. Here one should remember that it is electrons and holes that contribute additively to the current. 
The current is therefore determined by the charge carrier density, that is by the sum of the densities of electrons (n + ) 
and holes (1 — nr) — i.e. by n + + (1 — n~) = 1 + 2pj ) — rather than by the difference n + — (1 — n~) = — 1 + 2po. In the 
presence of a magnetic field one needs to involve po. ln the Hall component of the conductivity it is the difference of 
electrons and holes that contributes. 

The Boltzmann equation in matrix form is given by 

i[H, /] + d t p + -{vi, d Xi p} + eEid ki p - e zij eB z -{vi, d k] p} = J[p] (13) 



where the matrix- valued functional J is the collision integral. With the definition ( 12 ) the Boltzmann equation ( 13 ) 
can be written in the appealing form 

d t n + d x ■ j + e d k ■ (nE + j X B) = 2 Jo, 
2(a x by + d t s» + d x -j^ + ed k - (a»E + j" x B) = J„ . (14) 

Semiclassical kinetic equations deal as above with densities in phase space. Densities in real space are obtained by 
integrating the phase space densities over momentum, e.g. 

f d 2 k 

j(x,t) = J -—^j(x,k,t). (15) 

The Boltzmann equation is typically written in terms of the quasiparticle distribution f aa i(x,p,t) rather than in 
terms of the Wigner transformed density matrix p aiT i(x,p,t). The electron distribution p is the quantity in terms of 
which the current and densities are defined. The quasiparticles described by / on the other hand represent the free 
particles that satisfy the Pauli principle. Thanks to this the equilibrium state can for / be expressed simply in terms 
of the Fermi-Dirac distribution. This is of practical relevance in the analytical solving of the Boltzmann equation by 
linearizing around equilibrium. In this paper we neglect this difference. Thus, all expressions including p are assumed 
to apply for /. 



IV. DERIVATION OF COLLISION INTEGRALS 



The presence of, for example, electron-electron interactions, phonons or impurities, is captured in the collision integral 
J. It is assumed that one is in the kinetic regime, where the de Broglie wavelength Ab = 1/kp is much shorter than 
the scattering length £, which translates into £kp <C 1. We will only deal with averaged non-magnetic impurities. The 
averaging over of impurity positions in Vkk' = Y] n e~ l ^ k ~ k ^ Xn Ukk' restores the translational invariance.[69j 

In this section we will go through some of the possible approaches for the derivation of semiclassical kinetic equations 
including collision integrals from quantum theory. One group of approaches deals directly with the density matrix 
p(x\,X2,t) and starts from the von Neumann (quantum Liouville) equation id t p = [H,p]. The second type of 
approaches are Green's function techniques which have the Kadanoff-Baym or Keldysh equations for the double- 
time correlator G < (x\,ti,X2,t2) as their starting points. To recover a time-diagonal kinetic equation in the latter 
requires some approximation. This is the problem of Ansatz. We will see several different candidates and therefore 
several different collision integrals. One of them is identical to the collision integral derived with the density matrix 
approaches. 

A. Iterative solution of the von Neumann equation (quantum Liouville equation) 

The simplest derivation of a collision integral is probably one in which the von Neumann equation is iterated to second 
order in the interaction. Let H to t = Hq + V where V is an interaction switched on at a time to in the remote past. 
The von Neumann equation in the interaction picture is 

id t/ ? = \y\t),(?] (i6) 

withA 1 ^) = e lHot A(0)e- iHat = U f (t)A(Q)U (t). This is easily integrated to give 

p\t) = p\h)-i f e£\y\H),^H)\ (17) 
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which, when inserted back into (16), yields 



dtp 1 = -t[V l (t),p l (t a )]- / dt'iV 1 ®, [V(t'),p l (t')]} = 



= -i[V\t),p l {t Q )]- dt'{V I (t),{V(t'),p\t)}]- i / dt'iV^iV^t'UVit"),^")}] (18) 

J io J to J t 

after a first and second iteration, respectively. So far the equations are exact. The Born approximation allows us to 
get a closed equation for p at time t to second order in the interaction V . To this end we remove the last term in the 
second row. Alternatively, in the last term of the first row replace the full evolution with the free evolution, i.e. let 
p l (t') = p l (t), which has the appearance of a Markov approximation. Back in the Schrodinger picture the assumption 
of free evolution reads 

P (t>) = U (t',t)p(t)bll(t',t) = e -^°(*'.')p(i) e ^°(*'.') (19) 
and after reorganizing evolution operators one obtains the kinetic equation in the Schrodinger picture, 

d t p(t) - i[H ,p(t)] = -z[y, e - H «( t - t »)p(t ) e ^( t - t °)] - / ° dr[V, [e^Ve lH «\p(t)}}. (20) 



So far, this kinetic equation is locally time reversible (globally not, since we switched on the interaction). To capture 
the decoherence due to other processes (phonons etc.) we do not want the state to depend on correlations in the 
remote past. Therefore we include a factor e~ ,,T in the integral to impose this loss of memory. This introduces time 
irreversibility. This factor also regularizes the integral and allows us to send to — oo. In translating the evolution 
operators into Green's functions (see appendix [A]), the last term can be written as 

J[p] = -J ^[V,[G™VG° A ,p]}, (21) 

where we anticipated that the last term will become the collision integral J . 

Another source of irreversibility comes with the impurity averaging procedure, a coarse graining that also captures 
the decoherence due to phonons, for example. Then V in the first term becomes just a number ~ V(r = 0) and 
the commutator vanishes. |70j For the terms linear in ni mp one finds for example (summation over repeated indices 
implicit) 

(VG m VG 0A p) kk , — > n imp S(k -h + h- k 2 )u kkl G° k y klk2 G k K p k2k , = Z*G° k A p kk ,, (22) 
where we introduced the retarded self-energy = rii mv {uG^u) k . However, in a term like 

(VpG m VG 0A ) kk r — > n imp 6(k -kx + k 2 - k')u kklPklk2 G° k y k2k ,Gl A (23) 

the delta function seems to offer no simplification at all. 

At this point we can attain further simplification if we say that in the collision integral we are not interested in any 
contributions which have to do with non-diagonality in momentum. This actually amounts to saying that we are not 
interested in any gradient expansion corrections to the collision integral: 

(VpG m VG 0A ) kk — ► n imp u kk ,p k ,G k ?u k , k G° k A . (24) 

The collision integral can then be written as 



J[p(k,x,t)] = - W kk > J j^y(PkG° k K G° k A + GffiGl A p k - Pk'G™G° k A - G° k R G° k A p k > 



) (25) 



with the transition matrix W kk > — 27mj mp |w fcfc ' | 2 for spinless impurities. 

A quantum Liouville approach was also used by Culcer et al. [22 , although along different lines. The focus in 
that treatment was, from the start, only the diagonal part f k (not to be confused with our / for the quasiparticle 
distribution) of p kk i — f k 5 kk > + g kk > , closing the door to gradient expansion corrections in the interaction terms. The 
analogue of the iterative solution of the Liouville equation in this section is their decomposition into two coupled 
equations for f k and for the purely nondiagonal part g kk i, the integrated equation of the latter then being inserted 
into the former. Until here the approaches are equivalent. The difference comes with the Markov approximation. 
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Culcer et al. use f(t') — > f(t) as opposed to /*(£') — > With p(t') — > we find that the evolution operators 

cancel each other out in a different way so that at the end they sit around the entire inner commutator rather than 
only around the inner V, 



J[p] = - J ^[V,G™[V,p]G 0A ] 



(26) 



This is indeed the result in eq. (4b) in ref. 



We will see that the difference between (21) and (26) matters for 

Their 



the first quantum correction. We will also understand why it does not matter for the treatment of ref. [2 
recursive analysis taken to order (£k-p)° only requires the part that is insensitive to the differences between 
p6). 



21) and 



B. Nonequilibirum statistical operator approach 

The Nonequilibrium statistical operator approach (NSO) [47] is a full-fledged second quantized field theoretical for- 
malism and involves more conceptual and technical ingredients from statistical physics than the approach of the 
previous section. It was used in ref. 20J for the derivation of a pscudospin-cohcrcnt collision integral for graphenc. 
Like the Green's function formalism, the NSO formalism allows for Wick decompositions and therefore goes beyond 
the approach of the previous section for non-quadratic interactions. For quadratic interactions such as impurities the 
approaches should be equivalent. 

The starting point is again the Liouville equation 

id tP {t)-[H tot ,p{t)\ = - v (p(t)-p iel {t)), (27) 

where H to t in contrast to the previous section is the second quantized hamiltonian. r\ is a small number. The state 
variable p is the statistical operator. The density matrix is obtained as the expectation value of the relevant operators 
P m (to is a composite index). In particular, with P pcr p v = c^, a ,c va we have 



Pi' = (^,pv) ~ Tr(pP pff!pV ). (28) 



The novel ingredient on the right-hand side of (27) is included ad hoc to implement the irreversibility already 
on a fundamental level in the derivation of kinetic equations (thus well before the impurity averaging step in the 
previous section). This is done by introducing an auxiliary statistical operator /5 re i — the relevant distribution — that 
is the statistical operator with maximal entropy (and therefore minimal quantum mechanical correlations) among all 
operators with same expectation values for the relevant observables P m (to is a composite index) 

(P ro ) rel := Tr(p rcl P m ) = Tr(p rel P m ). (29) 



The relevant distribution p rc \ does not evolve according to a Liouville equation like ([27]) but is determined by the 
macroscopic observable (P m )t- It serves also as an initial condition p(t) — e~^ t-ti ) L p TO i(ti) where Lp := [H to t,p\- 
This evolution is however replaced by a coarse grained one, 

p(t) = — L- [ dhe-^-^p^ti, (30) 

t - to J t 

(eventually to — ¥ — oo) to reflect that a macroscopic system forgets the microscopic details of its initial state after 
some microscopic time. After some more steps for which we refer to the books by Zubarev et al. [47] one arrives at a 
general kinetic equation, here to second order in the interaction, 

^ + • ■ • = Jm(t) = jU(t) + ji 2 Ht) + 0(V 3 ) . (31) 
The first term of order V 1 is the mean-field contribution 

Jin\t) - ~ i([Pm,V])i el = -TV Mt)[P m ,V\) . (32) 
The second term is of order V 2 , 

S4n\t,t')' 



Ji 2) (t) := - f dtV^Tr L el (i') F,z4(M0[F,P m ]Z4(M0 + *E P ™^ 



\t' 



(33) 
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The first term is the familiar double commutator with the interaction. Notice that the average is taken with respect 
to the relevant statistical distribution, which allows for Wick decomposition (see [47]). For a two-body interaction 
one will find that the second half, the term with J^ 1 ' , does not contribute to the collision integral but cancels some 
anomalous terms from the first half. For impurities, J 1 ^ 1 vanishes after impurity averaging. Together with Born 
approximation (Markov approximation) 



p rel (t') = U {t',t)p^{t')ul{t',t)+O{V) 
which replaces the full evolution by the free evolution one arrives at [Uq, shorthand for U (t,t') 



(34) 



4? (*) [p] = - I dt' e "(*'-*>Tr (j4p Iel (t)U W, U\ [V, P m }U ] ) 

J —OO 

dt' e "(*'-*)Tr (p iel (t) [u VUl[V,P r , 



G VG , [V, p] 



(35) 



Notice that this is different from both (21 1 and (26). We will see that when the collision integrals for graphene is 
written out explicitly in spin components the result (21 1 and the result (|35|) will coalesce. 



We would like to mention that ref. [47] also treats the Cluster expansion method, which is advantageous compared to 
NSO for going to higher order in the interaction (deriving the full T-matrix), whereas NSO has an advantage for dense 
quantum systems (for deriving two-body collision integrals including the Pauli blocking factors). To second order in 
the electron-impurity interaction and to linear order in distribution functions both should apply equally. Nonetheless, 
the collision integral eq. (4.2.92) of ref. [27] is not equivalent to |2l| ) or ( [35] ) but to the result p6j). However, there 
is again a Markov approximation [see the section between eq. (4.2.24) and eq. (4.2.25)] where the full evolution is 
replaced by no evolution rather than by free evolution although one is in the Schrodinger picture. Correcting this 
point one obtains instead the result (21) with the Cluster expansion approach. 



C. Green's function approach 

In the Green's function approach one starts with general dynamic equations for the two-time correlator G < (ti,t2}- 
However, solving such equations is generally difficult and therefore some approximation that limits the equations to 
the the time diagonal t2 = t\ is desirable. This is also necessary if one wants to derive a Boltzmann type equation 
for pit\) — G < (t\,t\\ We start with the discussion of the Ansatz and then turn to the derivation of semiclassical 
kinetic equations from the Kadanoff-Baym equations. However, we will see that even for a given Ansatz one can 
derive different collision integrals. 



1. The problem of Ansatz 

The first proposed Ansatz was the Kadanoff-Baym (KB) Ansatz [3"T] 

G < (x,p,t,u>) = p(x,p,t)A(x,p,t,uj) (36) 

with A = «(G R — G A ) being the nonequilibrium spectral function. This is a slight nonequilibrium modification of 
the equilibrium result G < (fc,w) = /fd {<jS) A(k,uS) (the fluctuation-dissipation theorem) and therefore it is expected 
to be a good approximation close to equilibrium. For weak interactions one uses the quasiparticle approximation 
A 27r(efc — lu). 

For nonequilibrium beyond linear response the KB Ansatz fails. This was noted by Jauho and Wilkins 51 53 in 
Boltzmann treatments of transport in strong electric fields, where their results differed from those[Sl] derived with 
density-matrix methods. A similar discrepancy was observed in the linear conductivity when comparing with Kubo 
formula calculations |55j . Later Lipavsky et al. showed that the discrepancy could be cured with the modified Ansatz 
and coined the generalized Kadanoff-Baym Ansatz (GKBA) [38] (see also [30j I39j) 

G < (xi,ti,x 2 ,t 2 ) = i J d 2 x 3 (G R {xi,ti,x 3 ,t 2 )G < (xi,t2,X2,t2) - G < (xi,ti,X2,ti)G A (x 3 ,ti,X2,t 2 )) = 

= iG R (t 1 ,t 2 )p(t 2 ) - ip(t 1 )G A (t 1 ,t 2 )) (37) 
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(spatial variables suppressed in the second line) which reduces to the KB Ansatz in equilibirum. Ref. [55J showed that 
the right-hand side is the first term in an exact expansion, which makes it possible to address the range of validity of 
the Ansatz. The exact expression respects the causal structure of the Kadanoff-Baym or Kcldysh equations and also 
fulfills some other natural criteria. (See appendix [B]). Semiclassical gradient expansion and electric field modifications 
can now be treated in a more consistent way. 

The GKB Ansatz seems to be most common alternative in applications when the KB Ansatz is considered insuf- 
ficient. Interestingly, however, in general it does not give the same Boltzmann equation as the one derived with the 
mentioned density matrix approaches (Liouville equation approaches) . For the first quantum correction of graphene 
the difference matters. 

However, we do find an Ansatz for which the kinetic equation obtained with a density matrix approach is also 
obtained from a Green's function approach, namely if the GKBA is replaced by the anti-ordered version (AA for 
anti-ordered Ansatz) 

G < (t 1 ,t 2 ) = iG < (t 1 ,t 1 )G R (t 1 ,t 2 )-iG A (t 1 ,t 2 )G < (t 2 ,t 2 ) = t P (t 1 )G R (t 1 ,t 2 )-iG A (t 1 ,t 2 )p(t 2 ). (38) 

Although this Ansatz violates the causal retarded-lesser-advanced structure of KB equations and the Langreth-Wilkins 
rules [SB], it can be derived in a similar way as the GKBA (see appendix [b]) . The full result (including the omitted 
expansion terms) fulfills almost all the criteria required in ref. [55], in particular the causality requirement. The average 
of GKBA and the AA gives a third alternative, here named the symmetrized Kadanoff-Baym Ansatz (SKBA), 

G < (t 1 ,t 2 ) = ~(j>(ti)A(ti,t2)+A(ti,t2)p(t2)) ■ (39) 

This Ansatz to zeroth order in gradient expansion appears for example in ref. [25j . 

Considering the importance that the issue of Ansatz has for spinless electrons in nonequilibrium beyond linear 
response, we believe that the issue should be even more important for graphene calculations beyond linear response, 
at least when electron-hole coherent effects have to be taken into account. 



2. The problem of identifying the collision integral 

The generalized Kadanoff-Baym equation [571 EHj reads 

G< = G R £<G A + (1 + G R £ R )G 0< (1 + £ A G A ) (40) 

where all products are to be interpreted as convolution products in real space/time and in spin variables. The retarded 
and advance components are determined by the Dyson equations ((G ) -1 — £ R )G R = 1 and ((G ) -1 — S A )G A = 1. 
The self-energies are to first order Born approximation given by 

£< = n iiap VG<V, 

£ R ' A = n imp (V + VG m ' A V) -+ n imp VG on > A V , (41) 

where we will neglect the mean field terms ~ V 1 as we are not interested in shifts of the total energy in this paper. 

The term containing G 0< in (40 1 plays the role of boundary conditions and vanishes when acting with (G R ) _1 from 
the left or (G A ) _1 from the right, 

(G R )- 1 G< = £<G A , 

G<(G A r x = G R S<. (42) 

In particular, the difference gives the Kadanoff-Baym equation in differential form, which is a double-time precursor 
of the time-diagonal kinetic equations to be derived. For our discussion we want to write it in two ways. The first 
equation, to be called Gl, is 

[i<9 t -iJ,G<] = £ R G< - G<£ A + £<G A - G R E< . (43) 

It identifies all the self-energy terms of order V 2 with the collision integral. This is what we think should be done for 
a comparison with the Liouville equation based approaches of the previous sections, where all terms of order V 2 were 
identified with the collision integral. The second equation, to be called G2, is given by 

[id t - H -ReE R ,G<] - [£<,ReG R ] = i{ImS R , G<} - i{Y,<, ImG R } . (44) 



13 



It is a frequently encountered starting point of Boltzmann treatments that consider renormalizations and other 
quantum corrections. [30l ESJ [56l |57] Of the self-energy terms only those on the right-hand side are considered as the 
collision integral, those on the left-hand side, in contrast, as terms renormalizing the free drift. (The term Re£ R 
shifts for example the zero of energy. Thereby it shifts horizontally the conductivity as a function of gate voltage, 
e.g. shifts the minimum conductivity away from zero gate voltage. In this context, see experiments [121 113].) For 
spinless electrons the commutators on the left-hand side vanish if one stops at zeroth order in gradient expansion. 
In this case one obtains the same collision integral to zeroth order in gradient expansion as with (43). In general, 



however, and in particular the case of spin, the self-energy terms on the left-hand side contribute even to zeroth order. 
We hope to address such renormalization corrections in future work. In the present paper we are mainly interested 



G>). 



in the alternative structures that might be obtained for a collision integral from the right-hand side in ( 44 1 . 

Both equations (43) and (44) also hold for G > . Within the Keldysh formalism the same equations are derived for 
G K = i(G< 



Notice that the quantum Boltzmann equation[5l ], obtained by gradient expanding (44) to first order, is a semiclas- 
sical kinetic equation in the variables (a;,p,i,u;). Integrating the resulting equation over the frequency (independent 
energy) w gives a Boltzmann equation. |71j In this sense we solve Boltzmann equations in this paper, not quantum 
Boltzmann equations. 



3. Different collision integrals 



With two different ways of writing the Kadanoff-Baym equation and three different kinds of Ansatz there are possibly 
six new collision integrals. One obvious question is, which of them corresponds to the collision integrals of the previous 
section. The second and independent question is, which one is appropriate for the problem of quantum corrections 
to the conductivity in graphene. 

In this paper we believe we are able to present an answer to the first question. It seems clear that the pertinent 



collision integral is derived from Gl, that is (43). The question is what Ansatz to choose. Interestingly, it is not the 
GKB Ansatz, but the AA that returns the collision integral (21). The GKBA would give 

£ R G< - G<£ A + £<G A - G R £< = -iVG 0R V P G 0A - iG m pVG 0A V + iVG m P VG 0A + iG 0R V P G 0A V + ...= 



G<£ A + E<G A 



G R £< 



= -i[V,G m [V lP ]G t)A ] + ... (45) 

where, in each term like E R G < = VG 0R VG 0R p ~ VG QR V pG 0A , we neglected the parts that contain two retarded or 
two advanced Green's functions since such terms vanish when one integrates over the frequency to obtain the collision 
integral, 

J = - f^[V,G 0R [V,p]G 0A }. (46) 



We call this collision integral GlwGKBA (w for with). It is clearly different from (21 ). Interestingly it coincides with 
{261. 



For the collision integrals G2 derived from (44) note, for example, that with the GKBA we obtain 
1 



{£ R -£ A ,G<} 



i{S<,G R — G A } = - ^[V,G m iV lP ]G 0A ] - i[V, [G OA VG 0R 



,Pl 



(47) 



We denote this collision integral G2wGKBA. It takes the form of an average of (46) and (21 ), in the latter, however, 
with retarded and advanced Green's functions swapped. The other possible col 
next section. 



lision integrals are presented in the 



COMPARING COLLISION INTEGRALS FOR GRAPHENE 



The different possible general collision integrals discussed in last section will now be summarized. When writing 



J[p] = - J the integrands (. 



vN & 



NSO 



Gl 


GKBA 


Gl 


AA 


Gl 


SKBA 


G2 


GKBA 


G2 


AA 


G2 


SKBA 



of the various candidates are given by 
[V,G m [V,p}G 0A ] 



l[V,G m [V,p]G 0A } + 
l[V,G m [V,p}G aA ] + 
l[V,G 0A [V,p\G 0R ] + 
\ (G2wGKBA 
[G 0R VG 0A ,[V,p]] 



[V,[G 0R VG 0A ,p]] 
UV,[G 0R VG 0A ,p]] 
hv,[G 0A VG 0R ,p}} 
l[V,[G 0R VG 0A ,p}} 
f G2wAA) 



(48) 
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We do not have to solve the Boltzmann equation seven times since the first two cases are sufficient to deduce all cases 
except the NSO case, which, however, will turn out to coincide with the GlwAA/vN calculation. For example, the 
second term of G2wGKBA is similar to the GlwAA result with retarded and advanced Green's functions swapped. 
A closer inspection (see appendix [A]) reveals that this swapping has no effect on the delta function part of the collision 
integral, but changes the sign of the principal value part. Therefore the delta function part of G2wGKBA is given 
by the delta function part of (GlwGKBA +GlwAA)/2, whereas the principal value part is given by the principal 
value part of (GlwGKBA -GlwAA)/2. One can decompose the principal value part of GlwGKBA in two parts X 
and Y, where X is the part that is invariant when one compares GlwAA with GlwGKBA, whereas Y is the part 
that changes sign. Then one can work out the principal value parts J v = ±J VX ± J VY for all the above collision 
integrals, with the relative signs of J vx and J VY , respectively, determined by the scheme 





X 


Y 


Gl 


GKBA 


+ 


+ 


Gl 


AA 


+ 




Gl 


SKBA 


+ 





G2 
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G2 


SKBA 









(49) 



Note in particular that for the collision integral G2wSKBA the principal value terms vanish completely. 

The products in (48) are still general convolutions. At this point the spin structure is left intact but the space 



and time variables are Wigner transformed and the products are then gradient expanded in these variables. We now 
assume that gradient corrections in the interaction terms can be neglected. We also replace the Wigner transformed 
density matrix p with the quasiparticle distribution /. These are two approximations which might be incorrect when 
calculating quantum corrections, but give the framework within which we want to make a first step and compare with 



previous work. Furthermore, we assume non-magnetic impurities, u' 
impurity averaging the collision integrals become 



<T(T 
kk' 



5aa>u kk > with W kk > := 2im ilap \u kk/ 1 2 . After 



GlwGKBA J[f(k,x,t)} = - f k ,W kk , f ^(G^AfC^ + G™AfG k A ), 

vN J[f(k,x,t)} = -/ fc ,^/^(/ fe G fc R G^_/ fc/ G«?G fc A + h.c.), (50) 

NSO J[f(k,x,t)} = -f k , Wkk ,JJ^(AfG™Gl A + G™GltAf), 



to second order in the interaction and to zeroth order in gradient expansion. The shorthand notations J =: Ik' 
and A/ = f(k, x, t) — f(k', x, t) were introduced. The retarded Green's functions is the non-interacting one taken to 
lowest order in gradient expansion. For the spin-orbit coupled case it is of the form 

G °* R = Eth^V' ^ := \^ + °-°h) (51) 

s=± w - e k 2 

with uj + = uj + i0 + and with Q = uj — <f> being the gauge invariant frequency variable. In the Wigner representation 
one has X A = (A R )t. 



The three collision integrals in (50) would obviously be equivalent if the ingredients G 0R , G 0A and / commuted with 
each other, as is the case for spinless electrons. In the general non-commuting case, including the case of spin-orbit 
coupling, the collision integrals appear to be different. However, the different forms do not necessarily imply differing 
results. We will find that to lowest order in (£fcp) _1 they all reproduce the Drude conductivity obtained with Fermi's 
Golden rule. To higher order in quantum corrections, however, an agreement is not at all obvious. 

Notice that the collision integral for non-magentic impurities should generally satisfy the property J, J[f(k 1 x, t)] = 
0, expressing that in real space (a;, t) the collisions cannot act as a source or drain of particles of a given spin state. 
For all collision integrals except the NSO integral this is manifest since the collision integrals change sign under the 
renaming of dummy variables k ■<-» k' . In the case of the NSO result this is not manifest at this level but the explicit 
collision integral derived for graphene will turn out to have this property. 

For further comparison and for the explicit solution of the Boltzmann equation we will write the collision integrals 



(50) explicitly in terms of the components /o and /. To streamline the lengthy expressions some more shorthand 
notation is introduced, x means that the quantity x depends on primed variables such as k', s' etc, whereas x 
correspondingly depends on fe, s. For example S' = |(1 + <r ■ s'b k >). Also, Ax := x — x', for example Ae = — ef/ 
and A(sb) =sb- s'b'. 
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Inserting ( |5l| ) into ( 50 ) gives a collision integral J = J s + J v consisting of delta function terms and of principal 
value terms. The principal value part terms J v are given by 



GlwGKBA J*[f\ = f k ,W kk ,^j: ss ,V(£) 

vN&NSO J p [f] = S k ,w kk ,±Y, SS 'V(id 

The delta function terms are given by 



ss'bx b' 
2 

ss'bx b' 



(-A/ + <tA/ ) 
(+A/ + o-A/o) 



(7 



^ x A/ 
sSjV£ x A ^ 



(52) 



Jo"!/] = - fw kk ' ~$>(Ae 

ss' 

•/fe' z „„, 



1 + ss'b • b' 



A/ 



1 + ss'B 



with the matrix B(fc, k') acting on / given by 

GlwGKBA B 



vN 
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B 


rest 


B 



= +b(b') T 

= -i>(b') T 

= +b{b') T -&(b) T 

= b'(b) T 



b'(b) T 

b'(b) T 
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B' 


= B 


B' 


= B T 


B' 


= B 


B' 


= B T 



A/ 

A/ 



(53) 



(54) 



The "rest" stands for the GlwSKBA as well as all G2 collision integrals. As stated above, their delta function parts 
are just the sum of GlwGKBA and GlwAA/vN. Taking into account the momentum dependence of the basis vectors 



b' = bcosNA9 - cs'mNA6, 
c! = c cos NA9 + b sin NA9, 



(55) 



one derives 



GlwGKBA 


vN 


NSO 


rest 
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Bb = 
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bcosNA9 + 
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Bb = 


b' 


B'b' = b 


B'b' 


= b 
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= b 




B'b' 


= b 


Be = c' 


Bc = 


c' 


Be = 


c cos NA9 - 


bainNA9 


Bc = 





B'c' = c 


B'c' 


= c 


B'c' = 


= — c 




B'c' 


= 



(56) 



Note the particular simplicity of the approaches GlwSKBA and G2. 

For studies where a spin-coherent Boltzmann equation is linearized in a small b (6p ^ ^f) the terms ss'B do not 
appear. Also, it should be safe to neglect the principal value terms if one is only interested in the response without 
quantum corrections. In that case, the delta function part of the collision integral is the same in all formalisms. |27j 
However, for graphene (ep = &f) we need the full collision integral to calculate the quantum corrections and the 
approaches therefore differ. However, a crucial simplification comes through e = ±6 due to e = 0. With Ae = sb—s'b' 
we obtain (henceforth we write V (1/x) as 1/x) 



E ss ^(Ae) 

Ess' sS(Ae) = Y, S8 ' s'S(Ae) 
E ss ,ss>5(Ae) 

Ess' Ae = Ess' SS Ae 
Ess' S Ae 
Ess' S ' Ae 



2(S(Ab)+S(b + b')) 


2(5(Ab)-5{b + b')) 


2(s5 + ETE 7 ) 
z Ul b+b' I 



(57) 



and therefore half of the terms in ( 53 1 vanish, leaving 

J = -J k 'W w { 5(Ab) (cos 2 M A/o + a ■ A( 
, i 



1 + B 



/)) 



j( 6 + 6 ')( sin 2MA/ + o -.A(i^/)) 



^^M<T-(b + Kb')xAf + -^ a -(b-Kl/)xAf}. 

Notice that the equations for the charge fo and the spin / are completely decoupled for graphene. 



(58) 
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To obtain the principal value terms for GlwGKBA, let k = — and for GlwAA/vN/NSO let k = +. The principal 
value terms of the other cases can be constructed out of these. The term with 6(b + b') gives only a contribution from 
the point k — — k' , i.e. k = k' = which we will neglect [75] since we consider kp ^ 0. (This term assures formally a 
continuity in the conductivity when kp = + — > kp = 0.) Thus, the delta function part we consider is 



J S [f\ = - I W kk > a(A&)!(/-/' + B/-B'/') 

Jk' 



(59) 



For the approach GlwGKBA and for all G2 approaches we obtain due to (56) 

J & [f] = - J W kk ,6(Ab)A (j^bA - J W kk ,5{Ab) l -A (c/ e + zf z ) . (60) 

for the components / = bf^ + cfz + zf z . In next section we will see that the matrix (1 + B)/2 is responsible for the 
additional chirality- induced spin-overlap factor cos 2 (iVA0/2) occuring in the intraband transition rates involving only 
the probability densities ((6± |/|b±) =/o±/t of energy eigenstates. This is how the Drude conductivity found with 
Fermi's Golden rule builds in the suppression of backscattering in monolayer graphene (N — ±1). The "transition 
amplitudes" involving the off-diagonal components (6 ± |/|&t) — fz ± */e (the Zitterbewegung components) are 
more elusive and beyond the reach of Fermi's Golden rule. However, for the GlwGKBA and all G2 approaches, the 
result ( 60 ) shows that the scattering of the off-diagonal components becomes very simple since it contains no angle 
dependent chirality factors but only a factor i compared to ordinary spin independent scattering. 

VI. CONDUCTIVITY OF GRAPHENE WITH PRINCIPAL VALUE TERMS NEGLECTED 

In this section we calculate the electrical conductivity of graphene for non-magnetic impurities with the collision 
integrals J & given in the previous section. We neglect the principal value part J v '. We also assume low temperature 
so that — ep. For notational compactness we henceforth neglect the charge unit e, allowing it to appear only in 
final results. 

Before we start, for means of comparison, we will derive the Drude conductivity per valley by considering only 
one band (electrons /i > or holes fj, < 0), in which case we can apply Fermi's Golden rule. Let us consider 
electrons. For monolayer graphene the one-band Boltzmann equation linearized in the electric field (/ = / eq + 
with / oq = /fd {vpk — jj) for electrons) reads 



E-d k ,r = - [ S(vpAk)W kk ,cos 2 ^Af^ = -f^ [ 6(vpAk)W kk/ cos 2 ^-(l~cosAe), 

Jk' 1 Jk' 1 



(61) 



where the transition probability W kk > cos 2 ^ (contrast with W kk > for ordinary electrons) takes into account the 
chirality of the Dirac electrons, responsible for the suppression of back-scattering. The real space current is given by 

f I" kdk 

vf {E) = ^vpkT tr E k d k f FD (vpk - = Ea (62) 

k' J 8 J 27T 



where (reintroducing e and H to the right) 



|/i|r trF £k F e 2 

an ■= : = -; > 7rr"kp (63) 



with Ttrp := rt x (fcp), kp = \fj,\/vp and the mean free path I := uprtrp. In (62 1 we used the shorthand notation 



fe := 2?r J The result (63 ) is written in such a way that it also includes the case of holes (// < 0). To get the total 



Drude conductivity of graphene we multiply by a factor of four for the degeneracy in valley index and real spin. 



We turn to the coherent treatment of pseudo-spin and see that all approaches reproduce ( 63 ) to lowest order in 
'fcp) _1 but give different quantum corrections. With 

r°° k'(\k' 

K(k,9,e'):= / W kk ,5(b-b'), (64) 
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the compact notation E ■ d k — Efdje + Egj:dg, the decomposition / = f + cr ■ (6/g + c/g + zf z ) with deb — Nc and 



dec = —Nb and using the table (56 1 we find for GlwGKBA 



e ■ d k h 

N 

E ■ dkff, - Eg — fe, 
E-d k f £ + E § jf i) + 2bf z 
E ■ d k f z - 2bf & 



= X = 



7- 



= J = - 



J 



K cos 



2NM 



= / K [ cos 



2 

NA9 



? Af- b --smNA6(f £ + f e >] 



K sin 2 



NA9 



, (f £ + h , )--smNA9Af i> 



K^^Af. 



(65) 



The other collision integrals are obtained with the ingredients in replaced according to the following table 







j§ 




4? 


GlwG. 

vN 

NSO 

rest 


cos 2 ^A/ 
cos 2 ^Af 
cos 2 fA/ 
cos 2 ^A/ 


coB a ^A/ 6 -|sinJVAfl(/i + /«) 
cos 2iVAe A/g _ ismNMift-U) 
coB a ^A/ 6 -|sinJVAfl(^ + / e ) 
cos 2 ^A/ & - 1 sin NA6f e 


sin 2 ^ (/e + /.') + i sin NA9(f^ - / s ) 
cos 2 ^MA/ e +isiniVA^-/ 6 ) 
cos 2 ^MA/ e +isiniVA0(/£ + / 6 ) 
i/a - i cos JVA0/£ +i sin ^Afl(/£ - / £ ) 


sin 2 ^A/, 
cos 2 ^A/ 2 
cos 2 ^A/, 



(66) 



All terms including the trigonometric factor sinAiVS are in the assumed case of symmetric scattering K(A9) = 
K(—A9) actually the same in all approaches since only the part including /' can survive, whereas the part including 
/ vanishes trivially. Thus, the approaches differ only in the elements J7e[/e] and J z \] 7 \. These, however, will only 
enter the solution to order ~ (ikp)^ 1 . We will see that the first quantum correction to the conductivity will depend 
on J z [fz\. Since the iterative solution in Culcer et al. [22] was only taken to order {ik-p)° the choice of Markov 
approximation and in general the choice of formalism would not have mattered. 

For further comparison with ref. [22], notice that eq. (12b) for PiJ'fS'ii] (our Jc[f b \) and eq. (12c) for -PhJ^SjJ (our 
Jf\fc\) exactly match the GlwGKBA collision integral above. The comparison was, however, already established in 
the earlier stage (26). At the present stage, on the contrary, a comparison also with the NSO derived collision integral 
of Auslender et al. |20] is possible, see appendix [C] 

We proceed by linearizing the equations in the electric field with / = / cq + and by Fourier decomposing the 



components of fr = 2~2i 

Jo ±/g 
Since E ■ d k f ccl with 



mB frn , (f = 0, b, c, z). In equilibrium we have 
f FD (±v F k) = 6(/iT«Ffc) + 0((We F ) 2 ), 



fa 



req 



(67) 



-dkfrv (e±) 



Eg = -(e i9 £*-e- ie £) 



(68) 



(where £ := E x + iE y ) only contains n — ±1 Fourier components, we can right away conclude that frn — for 
n 7^ ±1. It is enough to study the equation for the n = 1 component because the n = —1 term is just the complex 
conjugate. We find for the n = 1 Fourier component of the linearized Boltzmann equation 



£* 
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(d k f c ^\ (i+ o 

1+ 

-il s 

I V o o 



. f/ e q 







+il s 

-2b 
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2b 
X K J 



( /of \ 

AE) 
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J el 

V /if } / 



(69) 



In terms of integrals defined below one has X K = X~ for GlwGKBA and 1 K = 1+ for GlwAA/ vN/NSO. For the 
GlwSKBA and all G2 approaches one has 1 K = (1+ +T~)/2. Likewise l x = 1+ for GlwGKBA, l x = 1+ for 
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GlwAA/vN/NSO and just the average for the rest. . 



for N = ±1 



1+ 


:= J gi Kcos 2 ^{l-e 


-iA0\ — „-l 
) — T tr 


= i/ e , 2f(l-cos2A0) 


1+ 


:= / ,^sin 2 ^(l + e- 




= 1+ 


x- 


:= / e ,^sin 2 ^(l- e - 




= Ifg, #(3-4cosA0 + cos2A0) 


il s 


f j£ sin 7VA<9 c ~iA8 

J 9' 2 


= | / e ,ifsiriiVA0sinA0 


= ±iZ+ 



(70) 



where we used that K(— A8) = K(A9). Notice that I s is odd in N whereas the other integrals I are even in N. 
The equation for fo is decoupled from the other components and is solved by /ox = ~^£* T trdkfo q with r tr 
. The other components are found by inverting the remaining 3x3 matrix, 




£* 



. JSJKjV -oq 



2|M| 



i2b(i s d k fl <1 



k J b 

N j-eq 

'b 



(71) 



^ fCq\ 



with the determinant \M\ = 4b 2 l + — (I S ) 2 I K . Adding up the two Fourier components n = ±1 one obtains 



(E) 



b 

■(E) 



E, 



(4b 2 +X A X" ) d k fp +X 3 X" f /_ 
4b 2 X++X+X A X"-(X 3 ) 2 X' t 

x«(z s a fc /r i +z+f / t eq ) 



fc (xr9 fe /g q + ^2 (<9fc/? q + £/j 



6( 4b 2 X++X+X A X»-(X=) 2 X K 

2b(X B fc /| q +X+f 
4b 2 X++X+X A X"-(X 3 ) 2 X'= 



NX" 



6 4b 2 



^0 2b 



1 rcq 
fc-'b 



(72) 



eq 



1 j-cq 



The far right shows the monolayer case TV = ±1. The result ( 72 ) is the solution to all orders in (ikp) -1 ■ The expansion 
can be carried out by expanding the denominator in 1/6. In the monolayer case N = ±1 the denominator simplifies as 
\M\ = Ab 2 I + due to X A = I + and I s = NX + , therefore no expansion is possible. Consequently, contributions beyond 
the leading order correction will be absent in the monolayer case. Our monolayer result has a structure similar to that 
of eqs. (13-14) in Trushin et al. [21j . however, up to the sign of the ];/~ q part. The sign difference, which originates 



from the sign of E^ f^ in the left-hand side of (65), will determine the sign of a 11 to be introduced below. 



The charge current j in momentum space, see eq. ([8|), is in the graphene case (eo = 0) given by 

j(k) = 2kf b d k b + 29f^ — > 2v F (kf B + N0f e ) 
with the monalayer case N = ±1 to the right. With J g kE^ = J g 6Eq = E/2 one obtains the current in real space 



(73) 



3 = 



! j{k) = (a 1 

Jk 



The conductivity is given by the contributions 

a I _ _ j kdk 



(4b 2 +l K l x )d k b+X"X B 



2tt 4b 2 X++X+X A X 
dk TV 



C dfc AT X " T 

J 2tt Jv 4b 2 X++X+X A X*-(X s ) 2 X 



pi, f°q 

(IS) 2 IK o fc / & 

b+X»X+^ 

■'b 



-%fkdk(^ + ^)d k fl« 



.«£ f Hfr 51 f cq 

2tt J QK 2b 2 •/ ft ' 



(74) 



(75) 



the 



where for the monolayer case to the far right it was used that d k b — | = vp for all k. Because of X S _ N = —X S N 
conductivity is invariant of the sign of N. 

As a part of a 1 we recognize the Drude contribution a — — J ^^^d k f^ q = ik-p/An. The contributions a 11 and 

5a 1 = a 1 — <Jo are quantum corrections, both of order (tkp)^ 1 and for \N\ < 2 containing also higher powers of 
1 L Notice that there is no contribution of the order (£kp)°. 

= one has from (67) that /^ q = — |©(fc — kp) and d k f? q = —\8{k — kp). The conductivity including the 
leading quantum correction ~ (Ikp)^ 1 [for general N obtained by a truncated expansion of the denominator in (75)] 
becomes 



Kp)' 

At T 
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where we write out the monolayer (N = ±1) result to the far right. With real spin and valley degeneracies included 
the conductivity of graphene is er grap hcnc — 4(fr J + cr 11 ). The leading correction depends on X K but not on X x . Thus 
the difference in X K — that is the one in J z [f z ] — is the one that is crucial for the discrepancy between approaches. 

The leading corrections are positive within all the approaches. When we later include principal value terms this 
will no longer be the case. 

For screened charged impurities W(k, A8) — 2nrii mp (2k sin -^7^ + fcTF)~ 2 the integral in a 11 is convergent. For point- 
like impurities K(k, A9) = jAjWo with W kk > = Wo = const oc ni mp the integral for a 11 has a logarithmic divergence 
in the monolayer case (since X oc k/v — k/v-p) unless an ultraviolet cut-off is introduced. Let us nonetheless make the 
observation that X + — K/4 and 1~ = 3K/4. Given a cut-off, the leading quantum correction is larger by a factor 
of 3 with the approach GlwGKBA (X K = I~) compared to the approaches GlwAA/vN/NSO (X K = Other 
approaches lie midway between these two results. 

For point-like impurities in the multilayer case \N\ > 2 all the approaches coincide because 

X + = -K=X + =X~ (77) 
hence 1+ = l x = X K =: X. Furthermore X s = 0. With b = ak N and X = \K — f^Wn the T = limit of the 



untruncated form ( 75 ) is easy to carry out. The result is (here with e and h reintroduced) 



a 11 - g2 N * .nrctnn |jV| { ' 

° — 2h 4(|2V| — 1) ALLld - 11 2ik F ■ 



All integrals converge without any ultraviolet cut-off. The N = 2 case of ( 78 1 was worked out in collaboration with 
M. Trushin et al. , see ref. ITOl. 



The treatment in this section would until equation ( 75 1 also hold for the density matrix p, provided that is 



independent of k and that /9~ q = = 0. The difference would come in the last step (76). If the difference between 



p? q and /? q is of order V r , r < 4 it could in principle change the value or the order of the leading quantum correction. 



For an ac field E(t) 



VII. AC RESPONSE 

e lut E one obtains with the Ansatz f k E \t) 



e %u)t Boltzmann equation 



Aujt 



(Wf> +i[H,fl B >}+E.d k f?-J[f l 



(£)- 
k 



from which the Drude result ( 63 ) is modified into 

£70 M 



£k F 



(79) 



(80) 



47r(T trF 1 + iu) 4tt(1 + iwr trF ) 

In the coherent treatment of spin of the previous section the term iwf^e 1 ^ 1 enters as a diagonal matrix 110X4, in 
equation (691, i.e. we obtain the ac result from the dc result for the n = ±1 Fourier components with the substitutions 



IUJ, 



IUJ, 



IUJ, 



X S ^X B , 



(81) 



Since the Boltzmann equation (79) with u> 7^ is no longer real, the n = — 1 component of the solution is no longer 
obtained by simply complex conjugating the n = +1 component of the solution. The correct n = —1 component 
is obtained by carrying out the substitution (81) after the complex conjugation. This the same as carrying out the 
substitution in the full dc solution where the two Fourier components have already been added up. With (81 ) the dc 
result (72 ) is modified into the ac result 



fP = - E kTM\ ( [ 4 & 2 + (Z X + ^){X* + iu;)] cV- cq +I S \X- + iu,] If, 



f- 
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-E^l" + iu;) (l s 9 fc /? q + [X+ + iu,] I/f q ) 



TP i f° q ^ 

^9 k J i 



E §J i ri 2b{l*d k f? + lX+ + iu]±f i 



1 t cq 



k J b 



OkJb 4b 2 -u> 2 > 

Tp 1 rcq 2b 
C 8i;Jj 4b 2 -LJ 2 ' 



(82) 



with the determinant \M\ also shifted according to (81 ). We complemented this with the pure sample limits (W kk > 
0) to the far right. 
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Ac terms are derived from dc terms by replacing a real quantity I by an imaginary quantity iu>. Thus the real 
frequency-dependent contributions in <r(u) step up or down in even powers of w/X, whereas an odd power would 
be needed to derive a (£kp)° correction in Recr(w) from the dc result CDrudo ~ £kp or its dc corrections ~ (ikp) -1 . 
According to this argument, there are no corrections (£kp)° to Kea(u). This result should be contrasted with the 
frequency-dependent corrections of order (£kp)° found by Culcer et al. (eq. (31) and eq. (27) in [2"2"]). 

For further details, see appendix |D| 

VIII. FIRST QUANTUM CORRECTION WITH PRINCIPAL VALUE TERMS INCLUDED 

In this section we include the principal value terms and recalculate the first quantum correction to the conductivity, 
this time in a recursive solution in the spirit of ref. |22j . The recursion is only taken to order (£kp) in the distribution 
function, but could in principle be reiterated to access terms of order (ttp) -1 and higher as discussed in appendix [e] 
However, the structure of the recursive setting makes it clearer why the (£k-p)° correction goes only into f z , whereas 
/gj fe. — and by consequence the current — get no contribution if principal value terms are neglected. When the latter 
are included the recursive setting clearly shows why ft and /g then can get corrections already of the order (^f) - 

We will start with giving the physical reason for taking principal value terms seriously in our semiclassical kinetic 
equations. General collision integrals derived within a quantum coherent approach typically contain principal value 
terms (a.k.a. reaction terms, off-shell terms, off-pole terms) alongside with the delta function terms (a.k.a elastic 
terms, on-shcll terms, pole terms). The delta functions convey the sharpness in energy of the idealized semiclassi- 
cal quasiparticles. The quasiparticles are the almost free particles that will distribute themselves according to the 
Fermi-Dirac distribution in equilibrium, whereas the electrons due to the interaction will be given by a distribution 
with fatter tails. |30| The principal value terms are a reminiscence of the quantum coherent nature of the underlying 
particles and captures the deviation from the classical point-like "billiard ball" picture conveyed by the fully semi- 
classical (i.e. quantum incoherent) Boltzmann equation. One such example is the principal value terms related to 
the quickly decaying coherences coming from the redressing of the quasiparticles within the interaction radius. [30] 
The corresponding decay time (the collision time, the spent within the interaction radius) is in the kinetic regime by 
assumption much shorter than the relaxation time (roughly Tt r ), wherefore the electron quickly recovers its asymp- 
totic quasiparticle nature after one collision on its way to the next. In the spectral function the off-pole part is 
the broad background around the quasiparticle peak. 58 For spinless electrons there are ways of separating out the 
off-pole part from the quantum kinetic equation, with the remains becoming the standard Boltzmann equation for 
the quasiparticles. [50] 

The electron-hole coherence (or spin-coherence), too, is a deviation from the fully semiclassical particle picture, in 
this case not because of interaction effects but because of the Zitterbewegung due to the spin-orbit coupling. It is 
therefore no surprise that a spin-orbit coupling contributes with its own principal value terms adding to those related 
to the quasiparticle redressing. However, this time we do not want to separate out the principal values in deriving 
a Boltzmann type equation since the Zitterbewegung is known to be inherent in the asymptotic free particle. When 
we want to derive a kinetic equation while keeping track of electron-hole coherent effects, we should therefore keep 
the corresponding principal value terms. To our knowledge this physical motivation has not been given before in the 
graphene context nor in related fields (e.g. spin Hall effect, anomalous Hall effect). The only reference we know of 
that treats principal value terms due to spin (although not spin-orbit coupling) is the Green's function derivation by 
Langreth and Wilkins|56] of a Boltzmann equation for electrons interacting with localized spins. There the principal 
value terms are found to give important renormalizations. 

The technical problem with including principal value terms is that the two momenta k and k' in the collision integral 
are no longer confined to sit on the same surface. This is in contrast with the previous situation which allowed us to 
plug out the Fourier coefficients f rn (k) from the integrals, of which the remains become matrix elements like I etc. 
(The problem of solving the Boltzmann equation to all orders in (^fcp) -1 then reduced to a matrix inversion.) With 
the principal value terms we have to confront difficult integro-differential equations. In ref. (20) an analytical solution 
to all orders in (£fcp) _1 was obtained for point-like impurities. The solution for screened charged impurities is still 
unknown and will be addressed here within the recursive scheme. It happens that we do not run into the problem of 
finding unknowns inside of the integrals ,7 P with principal values. Therefore, no simplifying assumption about the 
potential is needed. 



21 



From the principal value terms in (58) we find in the case GlwAA/vN/NSO 
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For a comparison with Auslender et al. [20], see appendix [C| The case GlwGKBA is obtained by simply exchanging 
the principal values 1/A6 and 1/(6 + 6'), as can be seen in 1§. Decomposing J p = J px +J PY for GlwGKBA into 
parts even and odd, respectively, under A6 ++ 6 + 6' we have 



pa 

b 



= o, 

= +4' 



7- 



PY 







1 




2tt 




fb' 






S- 


1 




2tt 




f&' ~ 





PY 



(/a - cos AA0/1 + sin NAO ft), J z 



PY 







/ 1 


1 


-J* 


2tt 


I b+b' 


Ab 






/ 1 


1 




2t 


I b+b' 


Ab 






/ 1 


1 


— Jfc' 


2tt 


\ b+b' 


Ab 



smNA0f z , 



(84) 



(cos N AO f e -f a ). 



The corresponding expressions for the other approaches considered i n th e paper are obtained by multiplying the J px 
and J PY components, respectively, with +1, —1 or according to (49). 

For the same reasons as before only the Fourier components n — ±1 of the nonequilibrium part f( E ^ can be nonzero. 
Therefore one only needs to consider 
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Both in (84) and in ( |85| ) have we left out all terms that vanish due to the assumed symmetry W(— AO) = W(A9) 
of the potential. We will see below that for the calculation of the first quantum correction one only needs to know 

Jj*[/s] and J[[f.], 'ha' Is 



(86) 



We can already at this stage extract some general conclusions for the multilayer case |JV| ^ 1 with point- like 
impurities. The trivial vanishing of the angular integral in (86) will imply below that principal v alue terms do not 
contribute to order (lfcp)°- Nor do they contribute to order (Ifcp) -1 (see appendix [e|) . Combining (491 with the fact 
the terms J PY in ( 85 ) vanish trivially we conclude that principal value terms do not contribute to any order for the 
G2 approaches. The integrals J px vanish non-trivially in the special case \N\ = 2, in which case the principal value 
vanish in all approaches (including the density matrix approach[10 ). To the orders to which principal value terms do 
not contribute, the quantum corrections are given by the treatment in section |VI| For the monolayer case principal 
value terms can contribute to the correction ~ (tkp)° in all approaches, with the trivial exception of G2wSKBA. 

The Boltzmann equation for (we henceforth drop the the superscript (E)) can be written as 



S[f]+J s [f}+J p [f], 



(87) 
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where D is the driving term with 
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S[f] is the spin-precession term 
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and the functionals J and J v are read off from (651, (66), (83 1 and the remarks below (83). A more informative 



way of writing the equation (87) is 
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For notational simplicity we now prefer to see the expansion of / in orders of (ik-p) 
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with oc W n oc (£kp)~ n . Here J*- -1 ) oc W~ x is the lowest order result that yields the Drude conductivity. Notice 
that the functionals J[f] increase the power in W by one whereas the action of S[f) is neutral in powers of W . 
Therefore, the two latter equations do not allow /g and f z to have a lowest order component W~ x , since S[f^^] 
would return a term of order W' 1 , which could not be matched by any of the other terms T> (~ W°) and J\f\ (~ W° 
and higher) . The absence of in the first equation (the diagonal components do not precess) is what allows only ft 

to have a term of order W- 1 . Solving V h = J b [f ( f 1] ] (~ W°), yields = -E k T tI d k f CCi . 

The components /(°) are found by solving the system 
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where f~ ^ is known. The two latter equations constitute a closed set, which allows us to first find f^ and /J . 

Only the known component /- ^ goes into the principal value part. The system of equation is solved as in the 
previous sections by Fourier decomposition and matrix inversion, 



iN£" j-oq 
2k Jf, 




(-1) 



+ 



(-i)i 



(93) 



Notice that if we discard the principal value terms J v we find that ff^ = 0, which using the first equation in (|92|) 



implies /i ' = 0. This is exactly what the solution ( pFlj ) tells us: there is no W° correction [i.e. (£kp)° correction] to 

the conductivity in (76); only f z obtains a contribution to order W°. The last line of (93) indeed corresponds to f z 
in d7lb. ' 



Including principal value terms yields a nonzero /~ ^ . It also gives a nonzero f^'' according to the first equation in 
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Nonzero in-plane components 
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result in a correction to the conductivity. In particular, this correction is of order (£&f) since the components in (95) 
are of order (£kp)°. 

A closer inspection shows that /£* is real and /£* is imaginary, as was the case in |7l| ). This implies [consult 



equations (|72j) and (|68|)] that /| 0) = Effi^ /£* and /i UJ = E § 2f£>/i£*. The first quantum correction 5a to the 
conductivity is therefore for arbitrary N given by 
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From ( 95 ) one obtains the first quantum correction as a sum of the contributions 
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where we introduced the notation Wkk'Ae := W kk >. It is easily seen that for point-like impurities the angular integral 
vanishes trivially for \N\ > 2. For further results on point-like impurities see appendix [E| 
At zero temperature the correction ( 98 ) can be written as 
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and in the monolayer case N = ±1 this can be simplified to (with d k b 
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where an ultraviolet cutoff k\ ^> fcp was introduced. Notice that only 5a x is ultraviolet divergent. 

The total quantum correction is 8a = 4 (±8a x ± 5a Y ) with relative prefactors (possibly zero) given by (49 ) for the 



different approaches. Notice that the leading quantum correction will not be positive in all approaches, in contrast 
to the situation in section |VT] For point-like impurities, in particular, we find 







Gf 


GKBA 


-log It 


8 


vN & NSO & Gf 


AA 




T 8 


Gf 


SKBA 


-log It 




G2 


GKBA 






G2 


AA 


T 8 




G2 


SKBA 








(102) 



The leading quantum correction is ultraviolet divergent for all the Gf approaches, including the density matrix 
approaches, whereas it is convergent for the G2 approaches. Furthermore, only the approach G2wAA gives a 
positive correction, namely 
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Ah 

In the Boltzmann regime Ik-p f this is a small positive shift to the much bigger Drude conductivity 4<tq = ^j-lk-p = 
\ for a constant potential V kk > = Vq. We mention that for screened charge impurities ultraviolet divergences 

" "■imp *q 

are absent. 

To obtain a contribution of order (ik-p) -1 , which should explain the initial onset of convexity in the conductivity 
as one approaches the Dirac regime, one iterates this recursive procedure, see appendix [EJ It is once again the case 
that one only needs to insert known distribution functions into the integrals containing principal values. Therefore an 
analytical solution is possible although increasingly cumbersome. By considering the expansion parameter (T + +iuj)/b 
the iterative procedure can be repeated in order to find corrections to the ac conductivity. The ac analogue is obtained 



with the substitutions (81) 



The result ( 103 1 applies to point-like impurities. Strictly speaking the assumption of negligible inter- valley scattering 
should break down, and it is questionable if the results can be used to discuss graphene experiments. This caveat 
does not apply neither for topological insulators with only one Dirac cone (see the end of the Introduction), (nor, 
of course, to numerical simulations of graphene including only one cone. Our main drive, however, is graphene with 
screened charged impurities and in particular monolayer graphene with the screening parameter q s := /ctf/^f ~ 3.2 
relevant for samples on silicon-oxide substrates. We will assume that this is already long-range enough for inter-valley 
scattering to be of secondary importance. In that case the two-valley results and one-valley results should be roughly 
the same, and we can discuss the former relying on results for the latter. We can now speculate that our leading 
quantum correction could be one of the contributions to the residual conductivity observed in the experiments of Chen 
et al. |f 3j . We plan to return to the quantitative analysis in future work. However, already here we can draw some 
qualitative conclusions based on the present section and recalling facts from sec. [TTJ The value of the correction can 
only depend on the dimensionless parameter q Sl which for monolayers is independent of hp and hence independent of 
the electron density. This leads to a rigid vertical shift of the Drude conductivity as a function of electron density 
as illustrated in fig. [f] The size of this shift depends only on natural constants and the dielectric constant present in 
&TF- Thus, the quantum correction could depend on the dielectric environment of the monolayer graphene sample. 

We have given a quantitative evaluation for the correction (ikp) in the limit q s — > oo because here this limit 
coincides with that of point-like impurities. (This was the situation for the Drude conductivity (see sec.|n|. However, 



there the absolute scaling with fcp in the scattering times was relevant, in contrast to the case of the correction (100).) 



We have also evaluated the shift ( f 00 ) in the opposite limit q s = < f of an unscreened Coulomb interaction. We find 
the corrections to be ultraviolet divergent and we find the sign of a Y to be the opposite. However, since q s ~ 3.2 > f , 
we expect the limit q s — > oo to be the more relevant limit. Therefore we expect also for the realistic value q s = 3.2 to 



encounter the case that only the one approach G2wAA gives a positive value and that value likely to be close to ( 103 1 . 
In case inter-valley scattering is negligible with q s = 3.2 this value should then also be relevant for the two-valley 
situation and thus for graphene experiments. 

Puddle formation due to charge inhomogeneities leads to a variation in the Fermi level, see e.g. ref. [50] . However, 
our shift should be insensitive at least to small variations as it is independent of the Fermi level. In the case where the 
impurities sit at a non-negligible average distance d from the graphene plane this introduces a second dimensionless 
parameter kpd, which depends on the density. In this scenario 8a becomes density dependent and the shift is no more 
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rigid. However, in the Fermi momentum range l~ x <C kp "C d -1 , which is the range where the Drude conductivity 
should be linear, the effect on the statements above should be negligible. Thus, the residual conductivity we make 
predictions for should be fitted only from the "strictly" linear part of the Boltzmann conductivity. 

Next we discuss the leading quantum correction (99) in multilayers. We have shown that for point-like impurities 
it vanishes trivially. Thus — referring to sec. [Tl] — we expect the correction to vanish close to the Dirac point. However, 
the interesting limit is the one far away from the Dirac point where ( 99 ) is all what remains of the studied quantum 
corrections. Here, it is more relevant to compare with the uncscreened Coulomb interaction. Irrespective of whether 
or not the correction is finite in this limit, it becomes doubtful that our analysis is still valid in this regime. The 
kinetic approach assumes that the collision time set by the range of the potential is short compared to the relaxation 
time T tr . For point-like impurities this is certainly the case whereas for an unscreened Coulomb potential it is highly 
questionable. For a long-range potential it certainly becomes important to take into account mean-field effects and 
their renormalizations of the free drift, which we have neglected. These issues are discussed in the context of strongly 
interacting spinless Fermi systems, see ref. |30] , At the moment we do not know how to generalize these issues to 
systems with non-trivial spin. Therefore, we cannot say much about whether there could be non-vanishing effects due 
to electron-hole coherence for multilayers with charged impurities. Since finite effects of electoron-hole coherence far 
away from the Dirac regime are very counter-intuitive we find it likely that in a proper treatment of the multilayer 
problem they would vanish at very high densitites. However, in monolayers we expect such effects to survive, as we 
come to next. 

Monolayer graphene stands out in many respects and comes with many surprises compared to multilayer graphene 
because of the linear dispersion and the unit winding number \N\ = 1. (See also refs. [2H [2Z] on why \N\ = 1 is 
special.) We saw already in sec.|H]that monoalyers are different to multilayers and 2DEGs when it comes to how the 
screening depends on the electron density. Therefore we should not straight away discard as unphyscial the finding 
of finite effects of electron-hole coherences far away from the Dirac regime, although we expect no such effects in 
general and in particular not in multilayers In both monolayers and multilayers the Fermi surface — and therefore 
the number of electrons contributing to a nonequilibrium response — grows linearly with kp. A Kubo formula for the 
conductivity (see e.g. eq. (2) in ref. [TO]) disfavors matrix elements between states with a big energy separation. 
The Zitterbewegung contribution from each electron would therefore be suppressed by the large energy denominators 
l/[e+(fc F ) - e-(k F )\ ~ k^ N . In the case of multilayers this suppression wins over the increasing density of states 
as k~F increases. However, in the monolayer case the two effects compensate each other, wherefore a finite effect of 
Zitterbewegung at large energy splitting is not inconceivable. 

One might worry about the electron-hole coherent effects being negligible compared to weak localization corrections. 
However, this is not necessarily the case, at least in idealized situations, as should be clear from recent numerics [10] . 
where the analytically found electron-hole coherent conductivity stays very close to the numerically exact value, with 
the small rather constant discrepancy probably due to weak localization. Nor should the electron-hole coherent shift in 
monolayers be negligible in the residual conductivity since we find it to be of the order of one quantum of conductance. 
Further, the different leading order quantum corrections can be cleanly separated and therefore treated independently. 
From the kinetic equation treatment of weak localization in ref. [49] we see that the weak localization correction takes 
only the Drude response part of the non-equilibrium Green's function as its input and not the full Green's function 
including contributions of higher order in (£fcp) _1 . Thus, like the Drude response, the weak localization correction 
should be independent of the choice of formalism. Therefore, we believe that the weak localization correction and our 
(ik F )° correction can be cleanly separated, and that the issue of formalism affects only the latter. 



IX. CONCLUSIONS AND OUTLOOK 



In this paper we investigated different derivations of semiclassical but spin-coherent Boltzmann equations in a case 
where differences could matter, namely in the electron-hole coherence originated quantum corrections to the Drude 
conductivity for 2d Dirac electrons, as encountered in in graphene or in the surface states of 3d topological insulators 
like Bii^Sba;, Bi2Te3, Sb2Te3 and Bi2Se3. With a few exceptions wc find these quantum corrections to be highly 
sensitive to the approach. We find the leading quantum correction in monolayer graphene to be particularly interesting 
as a litmus test, and suggest that a precise determination of this contribution from numerics or experiments might 
single out a unique approach. This sensitivity has motivated us to search for an Ansatz that provides the link 
between a quantum Liouville equation derivation and a Green's function derivation, that we find to differ with 
existing approaches. The simple structures of the derived collision integral in their most general form makes this 
search unambiguous. We have found the missing link, at least for the case of impurity interactions in the lowest Born 
approximation, and propose a novel Ansatz (AA) of a simple but counterintuitive form that to our knowledge has not 
been studied before. 

On a more technical level we pointed out that the fact that the pseudospin-orbit coupling is the dominant term in 
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the hamiltonian is essential for the differences to become important. The fact that spin-orbit coupling constitutes the 
entire kinetic part in the graphene case simplifies the collision integral considerably and makes an analytic solution 
possible. The analytic treatment becomes non-trivial due to the presence of principal value terms. We discussed the 
physical origin of these terms and explained why one also in our kinetic description should take them seriously. In 
addition, we showed how to deal with them for arbitrary scalar impurity potentials, at least for not too long-ranged 
potentials. We kept the winding number of the spin-orbit coupling general in order to address single layer graphene 
as well as multilayer graphene. 

We found that the first quantum correction depends both on the chosen formalism as well as on whether or not 
principal value terms are included. With principal value terms neglected the leading quantum correction is found to 
be of order (££:f) _1 . When they are included and do not vanish the leading quantum correction is of order (ikp) . 
An electron-hole coherence originated quantum correction ~ (^f) would be a counterintuitive result as it implies 
that electron-hole coherent effects could remain finite even far away from the Dirac regime. We discussed why such a 
result in the case of monolayers is not absurd, although surprising. In multilayers on the other hand we do not expect 
such a result, and indeed for point-like impurities the correction ~ (£kp)° vanishes trivially. For screened charged 
impurities one encounters the situation that the potential approaches the opposite limit — the unscreened Coulomb 
potential — when one increases the density. Since the kinetic analysis in this paper disregarded renormalization effects 
of the free drift, our analysis should break down when the quasiparticle spends a sizable fraction of its time within the 
interaction range of the impurities. Therefore, we do not attempt to evaluate the correction ~ (£kp) for multilayers 
with charged impurities in the high density limit where the screening is weak. 

We argued that the shift in monolayers due to electron-hole coherences should only depend on the dielectric constant 
through the dimensionless parameter q s . We also argued that the evaluated leading correction for monolayers with 
point-like impurities should be closely related to the one relevant to experiments since for monolayers on silicon-oxide 
substrates the screening parameter q s ss 3.2 is bigger than one and therefore more related to the limit q s — oo than 
to the limit q s — 0. Such a shift could be one of the contributions in the residual conductivity observed in recent 



experiments [13]. Our contribution 5a given in (100) depends crucially on the approach to deriving collision integrals. 
With a precise measurement of the residual conductivity and a precise knowledge of other contributions {e.g. weak 
(anti-)localization) that one would need to take into account, monolayer graphene would offer an unprecedented 
setting for experimentally singling out one among all the approaches restudied or introduced in the present paper. In 
a comparison with numerics one would of course have an even more controlled setting. Wc plan to address the actual 
quantitative analysis and comparison with experiments and numerics in future work. 

The observed residual conductivity is positive. If also the contribution from electron-hole coherences would be 
determined to be positive, this should pick the approach called G2wAA as the unique alternative. Interestingly, this 
contains the Ansatz AA that is one of our original contributions. In such a case, the electron-hole coherent effects of 
graphene would require a kinetic equation that not even in principle could have been derived with existing theory. 

The more technical work on the Ansatz and the introduction of the Generalized Kadanoff-Baym Ansatz (GKBA)[38 
was prompted by the study of high-field transport for spinless electrons, see ref. [39]. It would be interesting to study 
the consequences of our Ansatz on transport beyond linear response. The two Ansatzes certainly differ when electron- 
hole coherent effects are important, but it is not known to us if there is a difference for spinless electrons in strong 
electric fields. 

As already mentioned in the introduction we believe that all the presented approaches might still have to be refined, 
in particular for the case of weakly screened impurities, which is relevant for multilayers with charged impurities an at 
high fermi momentum. The elaborate literature on Boltzmann transport in spinless electrons offers two clear directions 
for improvement that we believe could also be important for the quantum corrections due to the pseudospin-orbit 
interaction. 

1. A proper accounting for of all terms that could contribute to first order in the electric field. For example taking 
into account that also the noninteracting response functions G 0R and G 0A are modified by the electric field 
and that the gradient expansion of the self-energy terms to first order includes electric field contributions in a 
gauge invariant formulation. In the context of spin-orbit interactions this was discussed recently in ref. |27) . 
For electron systems with a trivial spin index these issues have been discussed for more than two decades in the 
context of high-electric-field transport, see ref. [35] for a review. See also refs. [Tfl [57] . 

2. A proper extraction of the quasiparticle part / in the kinetic equation for p (or G < ) and a proper incorporation 
of renormalizations of the free drift. In the process one unveils the qualitative difference between the electron 
distribution function p and the quasiparticle distribution /. The difference appears as a wave-function renor- 
malization factor and as an extra term containing principal values. The latter term is related to the quickly 
decohering off-shell motion from the quasiparticle redressing within the interaction radius. For single-band 
electrons these issues have also been discussed for some two decades, for example in the context of Boltzmann 
treatments of Fermi systems with strong two-body interactions (for an extensive review we refer to ref. [30j , for 
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the context of impurities, see ref. [55] )• Here it becomes important to recognize that / eq = /fd 7^ p cq - Only 
by properly separating out the coherences related to the quasiparticle redressing can one in a controlled way 
extract a Boltzmann equation that one solves by linearizing around an equilibrium described by the Fermi-Dirac 
distribution /fd- 

It is an open question as to how one can generalize these issues to situations with interband coherences. However, 
already at the level of treatment given in the present paper we believe the discussed approach to be very promising 
for understanding electron-hole coherent effects in the conductivity in graphene, with the recent results in ref. [10] 
being one example of this promise. 
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Appendix A: Translating evolution operators into Green's functions 

With Hk = eofc + <x • bfc we notice that 
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Appendix B: Details on the anti-ordered Kadanoff-Baym Ansatz 

In ref. [5S] the correlator G < is divided into the auxiliary correlators 

G<\t x M) = 0(t 1 ^t 2 )G < (t 1 ,t 2 ), 

G <a (i!,i 2 ) = 6(t 2 -t 1 )G < (t 1 ,t 2 ). (Bl) 



By acting on G <r with (G R ) 1 from the left, using the generalized Kadanoff-Baym equation (42) and then acting on 
the result with G R from the left one arrives at 



G<\hM) = ^G R (t 1 ,t 2 )G<(t 2 ,t 2 )+ I' dt ( 2 dt'G? it X? t ,G< t2 +G? it X<G$ t2 , 

J £ 2 J —OO 

G <a (h,h) = -iG<(t 1 ,t 1 )G A (t 1 ,t 2 )+ I* dt [\t'G? lt2 X<Gt; t2 +G< t2 ^ t ,G$ t2 . (B2) 

We complemented this with the corresponding result for G <a with (G A ) _1 and G A instead acting from the right. The 
first terms to the right sum up to the GKBA. The integrals are correction terms that fulfill several natural criteria. 1) 
On the time-diagonal t\ = t 2 they vanish, making the GKBA exact. 2) No integrals stretch to t = +00, i.e. the result 
respects the causality of the Kadanoff-Baym equations. 3) One can derive the same equations for G > , i.e. particle-hole 
symmetry remains. 4) The spectral identity G < + G > = i(G R — G A ) is still satisfied. 

The solution can be used to determine G < iteratively to the desired precision. Thus, the GKBA can be seen as 
the first term in an expansion in the interaction strength. However, in ref. |38) it is noted that on top of that the 
arguments of the self-energies Stc run over disjoint intervals, which makes the integrals even smaller and relates it to 
the collision time To(<C r tr ), the small time the particles spends within the interaction radius. 

We now copy this treatment but act with the response functions fro m th e opposite sides (i.e. acting on G < with 
(G R ) _1 from the right using the generalized Kadanoff-Baym equation (42) and then acting on the result with G R 
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again from the right). This gives us instead 

G <r (t 1 ,t 2 ) = iG<(t u h)G R (t 1 ,t 2 ) + y tl dt'(G R E < + G<I] A -G <r S R ) tlt ,G^ 2 , 

G <a (t u t 2 ) = -iG A (t 1 ,t 2 )G < (t 2 ,t 2 ) + y t2 diGf lf (S R G < + E < G A -E A G <a )i t2 . (B3) 

The criteria l)-4) are still satisfied. Note in particular that the causality is respected. However, the result is a bit 
more complicated and this time the variables t and t' in the self-energies no longer run over disjoint time intervals. 
Our conclusion is that the expansion can still be seen as one in the interaction strength but no longer as one in the 
collision time. However, when comparing Boltzmann equations derived using Green's function techniques with one 
derived with a Liouville equation approach, it is sufficient to be consistent to the given order of the interaction. 



Appendix C: Comparison with the Boltzmann eq. of Auslender and Katsnelson 

The Boltzmann equation for the Dirac cone K was derived by Auslender and Katsnelson [20] with the NSO formalism. 
For comparison we are going to translate our Boltzmann equation into theirs. (Similar to us they do not consider the 
terms with S(b + b').) In ref. [2D]the Boltzmann equation is written for the quantities D := 2/ — 1, N := 2ft + 1 (not 
be confused with our winding number N) and g := f z — i/e- In terms of these we obtain 

E-d k D 

Es 

E -d h N + 2-Q-!mg 
k 

E ■ d k g - &N - 1) - i2bg 



Linearizing (e.g. D = D + 6D(E)) the equations in the electric field and comparing with (47), (51) and (52) in ref. |2"01 
we recognize all the terms. Apart from some minor differences in prefactors the main difference is the SD in equation 
(52) is in our case 5N. In our Boltzmann equation the electron density /o (i.e. D) decouples from the spin density 
/ (i.e. N and g). This decoupling is also found in ref. [22] . 



~ j , W *X cos 2 — AD, 



smA0 . 
i — - — N 
4 

'icos 2 ^Ag+ S -^N' i cos 2 ^fAg+ s -^N r 
^A6 + ir(b + b') 



(CI) 



Appendix D: On the ac response with principal value terms neglected 



In this appendix we concentrate on the monolayer case N = ±1. With (81) the determinant becomes 

|M| = 4b 2 (l + + iw) + (1 + + tuj) 2 (l K + iw) - (l s ) 2 (l* + iw) = Ab 2 (l + M2I LVf iIK ^ lU3) ) 

\ Ab z (l + + ioj) J 

The conductivity derived from ( 82 ) is given by 

fcdfc 1 



(Dl) 



a n (u) 



-v F 



-Vf 



2;t j¥| a fc ig q (46 2 + (X K + lW )(2I++ lW )) 

hc\h 1 



(D2) 



Under the assumption w,I <C bp = (here T = 0) the Fermi surface contribution a 1 can be expanded in and 
' ~ 4&2(z++^) ( 1 ~ t " (2 4b + ^|i++L) +m) ) t0 y ield the leading order correction 



with \M\ 



a\uS) « o a (u) 1 + 



2^(21+ +iu)(2£ : + iu) 
4| M | 2 (I++ lW ) 



(D3) 



29 



As already visible in (82), in the pure limit I — > 0, there are no corrections in a 1 to a 



tj \jA 



For 2 7^ 0, a 1 has, 

in contrast to the dc results of the last section, a correction of order (£kp)°. However, this correction only affects the 
inductive part Ima 1 . For u <C I + and for k = + we find that cr / (aj) 



a 



j— • i For the contribution cr 7/ 

the above expansion of | A/ 1 ~ 1 is allowed since 6 only sweeps the interval [\fx\, oo). However, the integral is cumbersome. 
We only give the result in the clean limit 1 = 0, 



' ,/fd(^-m)-/fd(-&-m) (t->o) . 1 2|/x| 

do - - = i m 

46 2 - lu 2 32tt 



ui 



4tt 7 46 2 - uj 2 32tt 2|/x| - w 32tt|/x| 

In contrast to the corresponding result for 1^0, this result does not diverge in the limit uj — > 0. 



(D4) 



Appendix E: Details on a second iteration when principal value terms are included 



In the second iteration the Boltzmann equation ( 90 ) reads 



f(i)i 



= 

o = S & [h , . c ,, 5 

= S,[/W] + J^[/i 0) ] 



Ji\fi°\fi 0) ] + jn/n 



(0), 



W 2 ) 
W 1 ) 
W 1 ) 



(El) 



The driving terms that were still present in the first iteration — see the equation ( 92 ) — are completely absent in the 



second and higher iterations. The components 



(0) 



/f and /<°>- 



-which determined the correction of order {tk-p)° in 



the conductivity — are the known input in the Boltzmann equation (|E1|) from which one then extracts the components 

^f) -1 in the conductivity. 



f~^\ /a and /J 1 ^. The latter components determine the correction of order 



We see that generally also the contribution of order (tk-p) 1 depends on the principal value terms. However, their 
contribution vanishes again trivially in the case of point-like impurities together with |7V| ^ 1, as it already happened 
in the first recursion with the consequence that f~^\ and the contribution ~ {(kp)° to the conductivity are all 

zero. For this case we immediately find /- 



terms and therefore consistent with (71 1. In contrast, /. 
The pertinent expression for GlwGKBA is given by 



(0) 



(i) 



[j£ [fz °']] will depend on principal value terms. 



see (69)] independent of principal value 



Jllfz 



1 

2^ 



W, 



k' 



kk' 



2 NAB 



Ab 



,2 NA9 



b + b' 



(E2) 



according to ( 49 1 we decompose ( E2 ) into 



derived with A6 ■<-» (6 + 6') from the result (83) for GlwAA/ vN/NSO. For the comparison with other approaches 



i 



6 + 6' 
1 



-iA8 H 



Ab 



co S NA6(f zl -e-* Af> f zl ) 



(E3) 



The contribution J vx does not vanish for \N\ ^ 1 whereas J VY does. From (49) we deduce that principal value 
terms in the approaches Gl (but not G2) can contribute to /j . However, this does not help as this contribution is 
removed in the last step, 



(E4) 



since J?[fz ] for \N\ ^ 1 vanishes trivially for the only nonzero Fourier components n = ±1. 

We have understood that principal value contributions to f z do not matter. Only if they appear in fg there can 
be a contribution to the current. Inspection of (83) shows that the only principal value terms that does not trivially 



remove the n = 1 Fourier component in the casepVj ^ 1 are the term J z f[fc] and the already considered term J^[f z ]. 



Thus J z [f^p] could contribute to fY' > . In a third recursion, i.e. for quantum corrections of the order (£k-p) 2 it is 



f(2) 
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therefore possible that principal value terms contribute to the conductivity even in the case of point-like impurities 
in the multilayer case \N\ ^ 1. The possibility of leaving out the principal value terms requires therefore further 
investigation. 
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